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A Guide to Tables on Punched Cards 


A list of tables on punched cards, compiled by W. J. EcKERtT, was pub- 
lished in the October 1945 issue of MTAC [v. 1, p. 433-436]. Since then, 
the number of laboratories interested in mathematical computations has 
greatly increased and as a consequence many new and important tables 
have become available. 

Computing techniques, too, have undergone radical changes with the 
advent of the electronic “‘C.P.C.” of the IBM Corporation and the new 
electronic multipliers which are accessible to laboratories of moderate size. 
Then there are the large, high-speed digital machines which employ punched 
cards; they are continually adding punched card tables of general interest. 

In a sense, the need for the more elementary punched card tables, such 
as square roots, reciprocals, and squares has decreased—at least for those 
laboratories that have electronic multipliers. In fact, several of the larger 
laboratories that formerly kept files of elementary tables on punched cards 
have discarded them. There are, however, a considerable number of smaller 
laboratories that still need the elementary functions, and even the labora- 
tories of moderate size would prefer a tested, available table of the Riemann 
Zeta function, say, to generating needed values with equipment currently 
at their command. In addition, “optimum interval” tables [MTAC, v. 1, 
p. 173-176] have found favor in some quarters. A number of tables of this 
character have been added since 1945. 

When, at the request of the editors of MTAC, the present authors under- 
took to bring the list of punched card tables up-to-date, it soon became ap- 
parent that the wealth of new material made it advisable to revise the 
guide completely. 

It was felt that classification of material according to subject matter 
rather than location would be most helpful to the laboratory that needed a 
certain table. In this type of classification it would have been difficult to 
include the file number which a given source attached to a table. Such 
information is therefore omitted. If a required table is adequately described, 
even a large and busy laboratory will probably find it quite easy to locate 
its file. An exception was made in the case of the UMT depository, that with 
time may acquire much heterogeneous material which might be hard to 
locate without a precise file number. 

Quite often the user of a table may want to know the original source of 
the keypunched values; whether such values were taken from well-known 
books, or whether the laboratory computed the table on punched card 
equipment from first principles. Moreover, one would like to know whether 
the possessor of a table checked it, or acquired it from another laboratory and 
did not get around to testing its accuracy. The authors tried to elicit this 
information. In the case of many tables of the higher mathematical functions, 
they knew from their own experience where such tables originated. For the 
elementary tables, the statements of the reporting laboratories were taken 
at face value, without exhaustive checks. Memory for such information 
becomes nebulous over the years, and many laboratories probably failed to 
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get credit for their own work; perhaps a few were credited for efforts they 
did not make. It is hoped that readers will report any such errors when they 
notice them. 

The tables reported are all “‘unclassified.’’ Some of the largest laboratories 
reported only a small fraction of their tables of general interest. Since such 
laboratories no doubt have good and sufficient reasons for limiting the library 
they wish to make available to others, the authors made no attempt to list 
tables which the source did not specifically mention. However, most tables 
of general interest are eventually published and reported in MTAC and 
private correspondence with the laboratories involved may bring fruitful 
results in such cases. 

The authors acknowledge gratefully the cooperation of the following 
members of laboratories who submitted lists of their tables: M. ABRAMOWITZ, 
H. H. Arken, P. Armer, W. D. BELL, J. BELZER, J. M. BOERMEESTER, E. C. 
Bower, C. C. BRAMBLE, S. R. BRINKLEY, Jr., B. F. CHEYDLEUR, G. M. 
CLEMENCE, R. P. Coates, C. F. Davis, W. J. ECKERT, B. FERBER, A. D. 
FRANKLIN, C. C. GoTLres, R. W. Hammina, C. Hastincs, Jr., R. HENSEL, 
F. H. HoLLanper, C. C. Hurp, J. P. Ketty, E. C. KENNepy, W. A. Kitts, 
3d, H. Krart, D. H. Lenmer, M. Lorkin, E. Lunpt, H. MCALLISTER, 
A. OpLer, H. PoLacwex, G. W. ReiTwiesner, C. V. Ruzex, J. ScHILT, 
V. SCHOMAKER, J. SHERMAN, C. T. Tar and F. M. Verzun. Many of them 
went to considerable trouble in answering correspondence with the authors 
and in checking the listing of the items reported by them. The difficult task 
of preparing the typed manuscript from a preliminary card index was per- 
formed by Mrs. KATHARINE WARREN. She not only checked the typed data 
with painstaking care, but also gave much attention to maintaining uni- 
formity in style and format. 

To the extent that this guide reflects accurately the tables available, 
credit is due to the wholehearted cooperation of all the laboratories involved. 
The authors of this guide alone bear the responsibility for any errors of their 
commission in reporting the data, and they will appreciate information 
about any inaccuracies discovered by readers. 

The meaning of ‘‘sources”’ and the symbols employed are explained below. 


SOURCES 


The source of the original tables from which card entries were key- 
punched is indicated wherever it is known. A list of numbered “sources” 
with the addresses at time of publication of this guide is given below. Some 
of these are laboratories which retain key-punched values. Others, such as 
(66) and (40), are books of tables or other texts to which reference is made 
in the guide. Some laboratories, like (64), no longer have the cards credited 
to them. In such cases, the place where the tables have been deposited is 
indicated alongside the name. Reference details about tables or texts were 
sometimes taken from (40) without further verification. 


MEANING OF SYMBOLS 


The symbol x = 0(.01)5(.1)8; 5D or 5S means as usual that the argument 
ranges from 0 to 5 at intervals of 0.01, from 5 to 8 at intervals of 0.1, and 
that the entries are given to 5 decimals (5D) or to 5 significant figures (5S). 
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The symbol [(41)-(9), (7) ] in the source column indicates that the en- 
tries are based on values in source (41), and that both (9) and (7) key- 
punched them independently. Similarly [(41), (9)-(9)] means that all 
key-punching was done by (9); that some entries were taken from (41) and 
others were computed by (9). The symbol [(9)] implies that computing 
and key-punching were done by (9). 

The laboratories that possess a given punched card table are listed in 
the column headed ‘Available at.’’ Some laboratories have stated that they 
have verified, mathematically or otherwise, the accuracy of the tables. This 
fact is indicated by the letter C following the source number. 

The symbol A*f (or merely A*) indicates the k* forward difference of 
the function; 6*f (or 6*) indicates the central difference of order 2k. Simi- 
larly &*f; and A*f; will sometimes be abbreviated to 57, A, respectively. 


CLASSIFICATION 


The classification followed is that adopted by FLETCHER, MmLErR, & 
ROSENHEAD in An Index of Mathematical Tables [source (40) 1. However, 
since the number of tables listed here is not very large, many sub-classifica- 
tions have been omitted. For example, the trigonometric functions are listed 
as follows: 


7.1 Tables with argument in radians 

7.2 Tables with argument in decimals of a degree 

7.3 Tables with argument in centesimal measure (grades or parts of the 
complete circle) 

7.4 Tables with argument in degrees, minutes, and seconds 


On the other hand, sub-headings such as 7.11, 7.12 have usually been 
omitted. 

When many similar tables are listed under one heading, the grouping is 
according to the number of decimals or significant figures given in the entries. 
For example: all tables given to 5 decimals or significant figures will be 
followed by those given to 6 decimals, and so forth. 

Some items seemed to belong to classes not specifically listed in (40). The 
following additional classes have been added: 


25. Random numbers and functions corresponding to random arguments 

26. Tables relating to aerodynamics. A wide variety of specialized tables, 
all of them available on punched cards, are listed in the voluminous 
unclassified reports issued by (54). Mention of the existence of such 
tables is made, but too much space would be required to include 
them in this guide. The reader interested in aerodynamic tables may 
consult the reports distributed by (54). Some of the punched cards 
are available at (6). 

27. Tables relating to number theory 

28. Tables relating to astronomy 

29. Actuarial tables. Here, too, no specific details are given. Complete 
details can be obtained from the Society of Actuaries or from the 
IBM Corp., 590 Madison Ave., New York, N. Y. 

30. Tables relating to map projections and geodesy 
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Some items belong equally well to more than one category. Thus 2" 
belongs under (2.3) as a table of powers and under (27) as a power of the 
prime 2. In this instance the table has been indexed in both categories. In 
other cases, where a listing is given in (40) the listing is usually included only 
under the category assigned to the item in (40), even though another choice 
might have been more desirable from a logical viewpoint. 

An effort was made to unify the entries as far as possible. On the other 
hand, two very similar tables may differ in minor details which might be 
important to the laboratory operator. For the most extensive table of a 
type is not necessarily the most desirable one; a compact table may often be 
preferred to a large table, given to a higher accuracy. For this reason, tables 
which differed sufficiently in format or other details are listed separately. 
This is especially true in the case of tables of trigonometric functions, and 
in the tables of powers. In this connection, it was not always possible to 
indicate in one item all the data available on one card. For example, some 
laboratories have the functions x!, 1/x, x*, and x* on one card. On account of 
classification difficulties, these are listed separately. 

In the case of optimum interval tables, the number of cards comprised in 
the compilation has been indicated, when this information was available to 
the authors. Unless otherwise indicated, such tables are interpolable linearly. 


SOURCES 


1 Bell Telephone Laboratories, Murray Hill Laboratory, Murray Hill, N. J. 
2 Wm. D. Bell, Telecomputing Corp., 133 East Santa Anita Ave., Bur- 
bank, Cal. 
3 California Inst. of Tech., Cooperative Wind Tunnel, Pasadena 4, Cal. 
4 Central Statistical Laboratory, Union Carbide and Carbon Corp., Oak 
Ridge, Tenn. 
5 Computation Laboratory, Harvard University, Cambridge 38, Mass. 
6 Consolidated Vultee Aircraft Corp., Ordnance Aerophysics Laboratory, 
Daingerfield, Texas 
7 Consolidated Vultee Aircraft Corp., San Diego, Cal. 
8 Dow Chemical Co., Western Division, P.O. Box 351, Pittsburg, Cal. 
9 E. C. Bower, Engineering Department, Douglas Aircraft Corp., Santa 
Monica, Cal. 
10 General Electric Co., Apparatus Dept., Schenectady 5, N. Y. 
11 Massachusetts Inst. of Tech., Cambridge, Mass. (Office of Statistical 
Services) 
12 North American Aviation Corp., att: Charles F. Davis, Tabulating Unit, 
International Airport, Los Angeles, Cal. 
13 Ohio State University, Cryogenic Laboratory, Columbus 10, Ohio 
14 Rand Corporation, 1500 Fourth St., Santa Monica, Cal. 
15 The Texas Company, Research Laboratories, Beacon, N. Y. 
16 University of Toronto, Computation Center, Toronto, Ontario, Canada 
17 Thos. J. Watson Astronomical Computing Bureau, Columbia Univ., 
New York, N. Y. 
18 Watson Scientific Computing Laboratory, Columbia Univ., New York, 
N. Y. 
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19 Central Computing Service, Univ. of Cal., Berkeley 4, Cal., att: P. L. 
Morton, Cory Hall 

20 California Institute of Technology, Gates & Crellin Lab. of Chemistry, 
Pasadena, Cal. 

21 Ballistic Research Laboratories, Aberdeen Proving Ground, Aberdeen, 
Md. 

22 U.S. Naval Observatory, Washington 25, D. C. 

23 National Bureau of Standards, Computation Lab., Washington 25, D. C. 

24 National Bureau of Standards, Inst. for Numerical Analysis, Los Angeles 
24, Cal. 

25 U.S. Dept. of Interior, Bureau of Mines, Pittsburgh 13, Pa. 

26 U.S. Naval Ordnance Laboratory, White Oak, Silver Spring 19, Md. 

27 U. S. Naval Proving Ground, Dept. of Computation and Ballistics, 
Dahlgren, Va. 

28 U.S. Army Map Service, Corps of Engineers, Washington 16, D. C. 

30 Society of Actuaries, % Service Bureau, IBM Corp., 590 Madison Ave., 
New York, N. Y. 

40 FLETCHER, MILLER & ROSENHEAD, An Index of Mathematical Tables. 
London, Scientific Computing Service, Ltd., 23 Bedford Square, 
W. C. 1, 1946 

41 Barlow's Tables of Squares, Cubes, Square Roots, Cube Roots and Recipro- 
cals. 4th ed., London and New York, 1944 

42 Joun Topp, Table of Arctangents of Rational Numbers. NBS Applied 
Math. Series, no. 11, March 1951 

43 Marchant Methods. Marchant Calculating Machine Co., New York, N. Y. 

44 Vega’s Logarithmic Tables of Numbers and Trigonometric Functaons. New 
York, D. Van Nostrand Co., 1912 

45 Mathematical Tables from Handbook of Chemistry and Physics. 9th ed., 
compiled by CHARLES D. HopGMan, Cleveland, Ohio, Chemical 
Rubber Co., 1948 

46 Epwin CHAPPELL, Five-Figure Mathematical Tables. 1915, reprinted by 
D. Van Nostrand, New York 

47 H. ANpovErR, Nouvelles Tables Trigonométriques Fondamentales, v. 1, 3, 
Paris, Hermann, 1915, 1918 

48 G. N. Watson, A Treatise on the Theory of Bessel Functions. 2nd ed., 
New York, Macmillan, 1944 

49 L. Scuwarz, Luftfahrtforschung. V. 20, no. 12, p. 341-372, 1944 

50 K. Hayasnt, Sieben-und mehrstellige Tafeln der Kreis- und Hyperbel- 
funktionen. Berlin, Springer, 1926 

51 F. Catiet, Tables Portatives de Logarithmes. Paris, F. Didot, 1795 
(Tirage, 1819). (sub-tabulated by E. C. Bower) 

52 J. J. Asrranp, Hiuilfstafeln zur leichten und genauen Auflésung des 
Kepler’ schen Problems. Leipzig, Engelmann, 1890 

53 J.C. P. Mitier, Tables of I,,(x), Kn(x), and Y,(x). The University Math. 
Laboratory, Cambridge, England. To be published for the Royal 
Society (the last of the British Assn. for the Adv. Sc. tables). Typed 
manuscript made available by author to (23). 

54 Applied Physics Laboratory, Johns Hopkins University, Silver Spring, 
Md. Handbook of Supersonic Aerodynamics, Jan., 1949. Reports can 
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be purchased from Supt. of Documents, Gov’t Printing Office, Wash- 
ington, D. C. Vols. 1 and 2 now available. 

55 SamMuEL HERRICK, Tables for Rocket and Comet Orbits. Punched cards 
available at (24). Functions tabulated: C,(U), S.(U), X.(U); C,(U), 
Si(U), X,(U), where 


(a) U = E —sin E; CU) = 1 — cos E; SU) = sin E; X.(U) 
= Efor0 < U < 3.15 

(b) U = sinh F—1; C,(U) = cosh F—1; S,(U) = sinh F; 
X,(U) = F for 0 < U < 300 


Especially useful for finding the functions corresponding to very small 
values of U. 
56 C. Cranz, Lehrbuch der Ballisttk. V. 1, Berlin, Springer, 1925 
57 Applied Mathematics Panel, National Defense Research Committee, 
AMP Report 24-1, Sept. 1943, submitted by MERRILL M. FLoop, 
Princeton Univ. 
58 Harotp T. Davis, Tables of the Higher Mathematical Functions. V. 2, 
Principia Press, Bloomington, Indiana, 1935 
59 R. L. ANDERSON & E. E. Houseman, Tables of Orthogonal Polynomial 
Values Extended to n = 104. Research Bull. 297, Agricultural Exp. 
Station, Statistical Section, lowa State College, Ames, Iowa, 1942 
60 D. H. LEHMER: (a) Tables of Ramanujan’s 7r(m), reported in MTAC, 
UMT 101, v. 4, p. 162, July, 1950 
(b) Table of the sums of fifth powers of the divisors of n, 
n = 1(1)5000, reported in MTAC, UMT 114, v. 5, 
p. 28, Jan., 1951 
61 National Bureau of Standards, Tables Relating to Mathieu Functions. 
New York, Columbia Univ. Press, 1951 
62 T.S. Keiiey, The Kelley Statistical Tables. New York, Macmillan, 1938 
63 F. E. Fow.e, Smithsonian Physical Tables. 8th rev. ed., Smithsonian 
Institution, Washington, D. C., 1933 
64 Dept. of Engineering Research, U.C.L.A. [tables available at (24) ] 
65 H. W. HoLtappeE., Tafeis van e*. Groningen, Holland, Noordhoff, 1938 
66 British Assoc. for the Advancement of Science, Mathematical Tables, 
London, England, Cambridge Univ. Press 
67 Admiralty Computing Service, H. M. Nautical Almanac Office, Herst- 
monceux Castle, Sussex, England 
68 J. PETERS, Seven-Place Values of Trigonometric Functions for every Thou- 
sandth of a Degree. New York, D. Van Nostrand, 1942 ed. 
69 AMELIA DE LELLA, Five-Place Table of Natural Trigonometric Functions 
to Hundredths of a Degree. New York, John Wiley and Sons, 1934 
70 U. S. Lake Survey, Military Grid Unit, Dept. of Engineers [cards 
at (23) ] 
71 ARNOLD N. Lowan & Jack LADERMAN, ‘Table of Fourier Coefficients,” 
Jn. Math. Physics, v. 22, 1943, p. 136-147 
72 H. J. Gray, R. Merwin & J. G. BRAINERD, Solution of the Mathieu 
Equation, Am. Inst. Elec. Engrs., v. 67, 1948, p. 429-441 
73 UMT Fire, (Unpublished Mathematical Tables) % D. H. Lehmer. 
(When writing for cards, mention UMT reference number.) 
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74 R. M. Rosrnson, (a) Stencils for the Solution of Systems of Linear 
Congruences Modulo 2. See MTAC, UMT 120, v. 5, p. 85. (b) Table of 
Integer Solutions of |y* — x*| < x. See MTAC, UMT 125, v. 5, p. 162. 

75 J. BAUSCHINGER & J. PETERS, Logarithmisch-trigonometrische Tafeln mit 
acht Dezimalstellen. V. 1, second ed., Leipzig, Engelmann, 1936 

76 Royal Observatory at Greenwich, Nautical Almanac and Astronomical 
Ephemeris, London, H. M. Stationery Office, London W. C. 2 

77 U.S. Navy, Hydrographic Office, Tables of Computed Altitude and Azi- 
muth, Publication No. 214, Washington, D. C., U. S. Government 
Printing Office, 1940 

78 The University Mathematical Laboratory, University of Cambridge, 
Cambridge, England 

79 C. E. Van Orstranp, Tables of the Exponential Function and of the 
Circular Sine and Cosine to Radian Argument. Nat'l Academy of 
Sciences, 5th memoir, v. 14, Washington, U. S. Government Printing 
Office, 1921 

80 H. BRANDENBURG, Siebenstellige trigonometrische Tafel . . . . Leipzig, 
Lorentz, 1931 

81 P. WiypDENEs, Five-place Tables. Groningen, Noordhoff 

82 H. M. Nautical Almanac Office (L. J. Comrie), Planetary Coordinates 
for the Years 1800 to 1940. London, H. M. Stationery Office, 1933 





Source Available at Description of Tables 
2. Powers, real and complex; simple ra- 
tional functions of powers 
(Also see 6.19) 
[(21)] (11)C, (16)C, 2.1 x?: x = 0(.0001)0.9999; exact 
(21)C, (24) On same card with x’, x*, x5 at (24) 
[(18)] (18)C x?: x = 0(.01)2.25; 2D; on same card 
with 2x?, x3, x4 
[—(3)] (16), (24)C x?, A: x = 1(1)1000; exact; on same 


(16), (24)C 


card with Vx, ¥i0x, 1/x 
x*?:x = 1(1)9999 A, 6; exact; on same 
card with 1/x 
Available at (17)C for x > 1000. 


[(41)-(27)] (27) x*: x = 1(1)10,000; exact; on same 
card with x’. 

[(41)-(18)] (18)C x?: x = 1(1)9999; 7D; on same card 
with 1/x, 1/x? 

[(41)-(21)] (21) x?: x = 1(1)12,500; exact 

(23)C x*: x = 10(10)10,000; 7S or 8S 
On same card with x*,k = — 1, —2, 

— 3, 3, 3, ;, ; 

[(18)] (18)C 2x?: 2 = 0(.01)2.25; 2D; on same card 
with x?, x*, x4 

[(21)] (16)C, 2.2 x*: x = 0(.0001)0.9999; exact; on 


(21), (24) 


same card with x?, x*, x® 
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Source 


[(18)J 
[(41)-(27)] 


[(21)] 
[(23) ] 
[(18) ] 


C(11)] 
[(41)-(27)] 
[(21)] 
[(16) ] 
[(18)] 
[(26) ] 
[(22) ] 
[(18)] 
[(18)] 


[(14)] 


[(41)-(21)] 
[—(3)] 
[(41)-(18)] 
[(41)-(27)] 
[(41)-(21)] 
[(41)-(26) ] 
[(23)] 
[(41)-(11)] 
[(23)] 
C(17)] 
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Available at 
(18)C 


(27) 

(21) 

(11)C, (23)C 

(18)C 2.3 
(11)C 

(27) 

(16)C, 

(21), (24) 

(16)C 

(18)C 

(26)C 

(11)C, (22)C 2.4 
(13), (18)C 


(10), (18)C 
(14)C 

(16) 

(21) 

(16), (24)C 
(18)C 

(27) 

(21) 


(23)C, (24) 


(23)C 
(11)C 
(16), (24)C 


(17)C 


Description of Tables 

x?: x = 0(.01)2.25; 2D; on same card 
with x?, 2x?, x* 

x®: x = 1(1)10,000; exact; on same 
card with x? 

x*: x = 1(1)12,500; exact 

x?: x = 10(10)10,000; 7S or 8S 

x*: x = 0(.01)2.25; 2D; on same card 
with x*, 2x?, x? 

x*: x = 1(1)2500; 8S 

x4: x = 0(1)10,000; exact 

x*, x5: x = 0(.0001)0.9999; exact 


x®, x7: x = 0(.0001)0.9999; exact 

x?: x = 0(.001)0.999; p = 2(1)6; 8D 

2"; m = 1(1)100; exact 

1/x, A: x = 1(.001)4(.01)9.99; 6D 

1/x: x = 0.1(opt.int.)3; 6S; 1896 
cards 

1/x: x = 0.629(opt.int.)16.15; 6S; 
1196 cards 

1/x: x = 1(opt.int.)10.00810; 5S or 
6S; 577 cards 

1/x: x = 2(.01)6.6; 7D; on same card 
with log, x 

1/x;A,B;A = 3(A; + Avs), B = 38 

x = 0.1(.0001)0.9999; 7S 

1/x, —A:x = 1(1)1000; 7D; on same 
card with Vx, ¥10x, x? 

1/x:x = 1(1)9999; 7D; on same card 
with 1/x*, x? 

1/x, A(1/x): x = 1(1)10,000 for 1/x; 
x = 1000(1)10,000 for A(1/x); 7S 


‘1/x: x = 1(1)12,500; 7S 


1/x, —A: x = 1000(1)2000(2)4000- 
(5)10,000; 7S; on same card with 
Vx, V10x 

1/x: x = 10(10)10,000; 7S or 8S 

1/x: x = 1(.1)20 and x = 1(1)2500; 
8D 


1/x: x = 1(1)9999; A, &; 10D; on 
same card with x? 


1/x: x = 1000(1)9999; 10D; on same 
card with x? 
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[(23 
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[(43 
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[(18 
[(21 
[(41 
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[(41 
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Source 


[(26) ] 
[(41)-(18) ] 


[(23)] 
[—(3)] 
[(25)] 


[(43)-(13) ] 
[(18)] 


[(18) ] 


[(18)] 
[(18)] 
[(21)] 
[(41)-(21)] 
[(26)] 


[(26)] 


[(26) ] 
[(23)] 


[(2)] 


[(41)-(27)] 


[(41)-(21) ] 


C(11)] 
[(41)-(27)] 


[(23)] 
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Available at 


(26) 

(18)C 
(23)C 

(16), (24)C 
(25) 

(14)C 
(13)C, (18) 
(18)C 


(18)C 
(24) 


(18)C 
(18)C 
(23)C, (24) 
(21) 


(23), (26) 


(26) 


(23)C, (26) 
(23)C 


(24) 


(27) 
(24) 


(21)C 


(11)C 
(27) 


(24)C 


2.5 


2.6 


Description of Tables 


1/x, A: x = 1000(1)19,999; 10D 
1/x*:x = 1(1)9999; 7D; on same card 
with 1/x, x? 
x, x-*: x = 10(10)10,000; 7S or 8S 
Vx, ¥10x, A: x = 1(1)1000; 4D 
Vx:x = 1(opt.int.)100; 5S; 401 cards 
Vx:x = 1(opt.int.)100;5D; 1082 cards 
Vxix = 1(opt.int.)100; 6D; 332 cards 
Vx: x = O(opt.int.)0.9950; 6D; on 
same card with ¥i — x; 583 cards 
Vx:x = 0.5(opt.int.)2; 6D; 246 cards 
Vx: x = 1(opt.int.)99.430; 6D; 181 





cards 

Vx: x= 1(opt.int.)100; 6D; 1793 
cards 

Vx: 1(opt.int.).0001; to about 


xs= 

OD; 242 cards 

Vx, ¥10x, A:x = 1000(1)2000(2)4000- 
(5)10,000; 6D 

Vx, V¥10x: x = 1(1)12,500; 7D for 
x < 1000, 6D for x > 1000 

Vx, A, 8, 3(Ao+A4): x = 101- 

(1)10,000; 6D 

Range extended by (23) to include 
x = 1(1)100. 

x, Vi0x, A: x = 1000(1)1995; 2000- 

(2)4000(5)9995; 6D 

, A: x = 10,000(10)100,000; 6D 

, ¥10x: x = 1000(1)10,000; 8D in 
Vx, 7D in Vi0x 

Vx, A, 1/Vx, A: x = 1(1)9999; 6D in 
Vx, 8D in 1/Vx 

[last place inaccurate ]; on same card 
with 1/x, A(1/x) 

1/Vx: x = 1(1)1000; 7D 

1/Vx: x = 1.000(opt.int.)99.180; 7D; 
902 cards 

Vx, A, B;A = $(A;+ A), B = $8?: 

x = 0.0001(.0001)1(.001)9.999; 8D 

Vx: x = 0(1)1200; 8D 

Vx, Vi0x, A; x = 1,000(1)10,000; 
8S-9S; also Vx, x = 1(1)1000 

Vx, A, &:x = 1(1)1000; 10S 


— 


Ah 


R 
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Source 


[(23)] 


[(26) J 


[(18)] 


[(26) ] 


[(23)] 
[(26) J 
[(82)-(17)] 
[(25)] 
[(26) 
[(41)-(21) ] 
[(41)-(21)] 


[(26)] 
[(41)-(27)] 


[(23)-(24) ] 
[(23)-(24) ] 
[(23)-(24) ] 


[(23)-(24) ] 
[(18)] 
[(67)-(24) ] 


[(13)] 
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Available at 
(24)C 


(26) 
(18)C 


(26) 


(23)C 
(23)C 
(17), (18) 
(25) 

(23)C, (24) 


(21) 
(21), (24) 


(23)C 
(27) 


(16), (18), 
(23), (24)C 


(18), (23), 
(24)C 


(24)C 


(24)C 
(18)C 
(24)C 


(13)C 


2.63 


2.7 


2.8 


2.81 


2.9 


Description of Tables 


Vx, V¥10x; A, 8: x = 1000(1)10,000; 
10S 
Tables for reducing interferometer 
data 
f(x, 9) = (N(x + 1)? — »* — V(x? — y*)/ 
(2x + 1) 
x,y = 1(1)100; y < x;5D 
Table of f = Vi — x:x = O(opt.int.)- 
0.9950; 6D; 583 cards 
On same card with Vx 
(1 — x*)-?: x = 0(.005)0.1(.001)- 
0.808(.0005)0.9; 4D 
to 6D 
x = 0.9(.0001)0.96- 
(.00001)0.99499; 5D 
to 7D 
(1 — x*)-+: x = O(various)1; 8D 
x: x = 10(10)10,000; 7S or 8S 
x—!: x = 2(.001)7.5(.01)20; 8D 
xt: x = 1(opt.int.)1000; 5S; 540 cards 
xt, A: x = 1,000(1)2,000(2)4,000(5)- 
9,999; 6D 
At (26)C with slight modifications. 
xt, A, B; A - 3(A; + Ai-1), B = $8; 
x = 1(.001)9.999; 7D 
xt: x = 1(1)12,500; 7D for x < 1000; 
6D for x > 1000 
xt: x = 10(10)10,000; 7S or 8S 
xt: x = 1(1)10,000; also A for x > 
1000; 8S 
xt, (10x)#, (100x)#: x = 0(.01)10; 9D 
or 10D 





x4, (10x)#, (100x)#: x = 0(.01)10; 9D 


or 10D 
xt, (10x)#, (100x)#: x = 0(.01)10; 15D 
in x; 14D in (10x)? and (100x)! 
x?, (10x)!, (100x)?: x = 0(.01)10; 15D 


xV8:x% = 0.1(opt.int.)1; 6D; 284 cards 

xt, xt, xt, xt: x = 1(.01)10(.1)100 
(1)1000(10) 10,000; 10S 

10-"°10*, a = 0, 1,-2: x = — 6.99- 
(.01)0.99; 7D 

Powers of complex numbers 








[(16 
[(24 


[(10 


[(41 
[(21 


[(26 


[(24 


[(18 
[(23 
[(45 
[(18 


[(13 
[(23 
[(44 


[(44 
[(64 


[(64 


[(18 
[(18 








Source 


[(16)] 
[(24)] 


[(10)] 


[(41)-(27)] 
[(21)] 


[(26) ] 


[(24)] 


[(18)] 
[(23)] 
[(45)-(4)] 
[(18)] 


[(13)] 
[(23)] 
[(44)-(27) ] 
[(44)-(21) ] 
[(64) ] 
[(64)] 


[(18)] 
[(18)] 
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Available at 
(16)C 


(24)C 


(10) 


(27), (28) 
(21) 


(26) 


(24)C 


(18)C 
(23)C 

(4) 

(10), (18)C 


(13)C, (18) 


(23)C 
(27) 

(21) 

(16), (24) 


(24) 


(26) 


(10), (16), 
(18)C 


2.96 
3.11 


3.3 


4.4 


6.1 


Description of Tables 


Vw? +i = a + ib: w = 0(.01)17; 8D 
ina; 7Dinb 

(x + ty)" = tn + in: 
x = — 20(2)(—10)(1)0; 
y = 0(1)20(10)400; 
n = 0(1)10; exact 

1/e:2 = x + iy;x,y = — 4(.02)4;5D 

Factorials 

x!:x% = 0(1)100; 8S 

x! and 1/x!:x = 1(1)1000; 62S; Com- 
puted on the ENIAC in connection 
with determination of e 

Binomial coefficients 

n = 1(1)100; exact 

Modified Bernoulli polynomials 


b.(x) = > n~* sin [nx + 4(k — 1)x], 


n=1 
k>1 
= [(—1)*5(24)*Bi(4x/m) ]/k!, 
where B, is the Bernoulli 
polynomial : 
x = ry/36, y = 0(1)36; 
k = 1(1)11;17D 
Common logarithms and antilog- 
arithms 


logio x: x = 1(1)99; 3D 

logio x: x = 1(1)13,000(10)23,000; 5D 

logio x: x = 1(.001)9.999; 5D 

logio x: x = 1(opt.int.)10(10*); 
4 = —1(1)5; 5D; 1463 cards 

logiox: x = 1(opt.int.)10; 6D; 568 
cards 

logio x, A: x = 10(10)10,000; 7S or 8S 

logio x: x = 1000(1)10,000; 7D 

logio x: x = 1(.0001)9.9999; 7D 

logio x:x = 1000(opt.int.) 10,000; 7D; 
586 cards 


logio x: x = O(opt.int.)0.962400; for 
quadratic interpolation; 7D; 119 
cards 

logio x: x = 1(opt.int.)9.9958; 8D 


logio x: x = 1(opt.int.)10; 8D; 5640 
cards 
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Source 


[(64) ] 
C(75)-(2)] 


[(75)-(2)] 
[(18)] 


C(8)] 


C(8)] 


[(46)-(21) ] 
[(44)-(22) ] 


[(23)-(14) ] 
[(64) ] 


[(S1)-(18) ] 
[(23)] 


C(18)] 
C(18)] 


[(23)-(11) ] 
C(23)-(11)] 
[(26) ] 
[—(3)] 


[(23)-(18) ] 
[(23)-(18) ] 
[(66), 
(24)-(24) ] 
[(23)-(12) J 
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Available at 
(16), (24) 


(2)C, (18)C 


(2)C, (24)C 


(18)C 


(8) 


(8) 


(21) 
(22)C 


(14)C 
(24) 


(18)C 
(23)C, (24) 


(4) 
(18)C 
(18)C 


(11)C 
(11)C 
(26) 

(16), (23), 
(24)C 

(4) 


(18)C 
(18)C 
(24)C 


(12)C, (24) 


6.26 


6.4 


7.0 


7.1 


7.12 


Description of Tables 


legio x: x = O(opt.int.)0.998284; 8D; 
707 cards 

logio x, A: x = 1(.0001)10; 8D. (Pre- 
pared at Vega Aircraft Corp. by 
Ray CRAWFORD, in collaboration 
W. D. BELL.) 

logie x, A, 1/A: x = 1(.001)9.999; 8D 

Radix table 

f = 10**; m-10*, m = 1(1)9; 
n = 1(1)4; to about 7S; 9 cards 





logio (: — *): x = 0.02(.01)0.99; 


y = 0(.005)0.05(.01)0.2; 5D 

Antilogarithms 

107: x = O(opt.int.)0.999; 5D; 407 
cards 

107: x = 0(.0001)1; 5D 

107: x = 0(.001)0.999; 8D; auxiliary 
functions for interpolation 

107, A: x = 0(.00001)1; 8D 

107: x = O(opt.int.) —5.3467875; 9D; 
707 cards 

10-*: x = 0(.025)1; 10D 

107: x = 0(.00001)1; 10D 

Trigonometric functions. Also see 24.2 
and 28. 

sin x, cos x: x = 0(.001)6.283; 5D 

sin x: x = 0(.001)1.571; 7D 

sinx: x = 0°(opt.int.)90°; 8D; 176 
cards 

sin x, cos x: x = 0(.1)20; 8D 

sin x, cos x: x = 0(.001)1.599; 8D 


- sin x: x = 0(.01)99.99, A; 8D 


sin x, cos x; A, &: x = 0(.01)27; 8D 


sin x, cos x, tan x: x = 0.0001(.0001)- 
1.5709; 9D for sin x, cos x; 8S (or 
less) for tan x 


sin x: x = 0(.001)1.999; 9D 
sin x: x = 0(.001)1.571; 9D 
sin x, cos x, A, &:x = 0(.0001)1; 11D 


tan x, cot x:x = 0(.0001)2; mostly 8S 
on same card with tanh x, coth x 








c(i 


[(6 
[(6 


C(1 
[(6 


[(6 
(22 


[(6 
[(4 
[(6 
[(2 
[(2 
C(2 
C(2 
[(2 


L(2 


C(s 


c(é 
C(t 
C(s 





)- 
r 





Source 
[(18)] 


[—(3)] 
[(68)-(27)] 
[(69)-(21) 


[(11)] 
[(68)-(11) ] 


[(68)-(17), 
(22) ] 


[(68)—(21) ] 
[(47)-(22)] 
[(68)-(26) ] 
[(23)] 
[(23)-(24) ] 
[(23)] 
[(23)] 
[(23)] 


[(27)] 


[(S1)-(14)] 


C(81)-(1)] 
[(18)] 
[(S1)-(14) ] 


A GUIDE TO TABLES ON PUNCHED CARDS 197 


Available at 


(10), (13), 
(18)C 


(24)C 
(27) 


(21) 


(11)C 
(11)C, (24) 


(10), (13), 
(16), (17)C, 
(18), (22)C, 
(23)C, (24)C 


(21) 

(22)C 

(26) 

(23)C 

(24) 

(24)C 

(23)C 

(23)C, (24)C 


(27) 


(14)C 


(1), (4) 
(18)C 


(14)C, (24) 


7.2 


7.3 


Description of Tables 
y = tan f: f = 89°.942(.001)90°; 3D 
Note: This is used by (9) as a criti- 
cal table for f = arctanx in the 
region close to f = 90°. 
sin x, cos x, tan x, cot x, A: 
x = 0(°1)360°; 5D 
sin x, cos x, tan x, cot x: 
x = 0(°001)90°; 4S — 7S 
sin x, cos x, tanx: x = 0(°01)44°99; 
5D 
tan x: x = 0(°1)90° 
sinx, cosx, sinkx, &; k = 1(1)9; 
x = 0(°1)360°; 7D 
sin x, cos x, tan x, A: x = 0(°01)90°; 
7S or 7D 
[range extended to 180° at (18); 
argument also given in deg., min. 
and sec. at (22) ] 
sin x: x = 0(°001)89°999; 7D 
sin x, cos x, A: x = 0(°01)90°; 8D 
sin x, tan x, A: x = 0(°01)89°99; 8D 
sec x, csc x: x = 0(°01)90°; 9S 
sin x, cos x: x = 0(1°)360°; 15D 
sin x, A: x = 0°(°01)90°; 15D 
sin x, cos x, &: x = 0(°01)90°; 15D 
sin x, cosx, cscx, secx:x = 0(°01)45°; 
15D in sin x and cos x; 8D in csc x 
and sec x 
sin x, cos x:x = 0(°01)90°; 17S. Com- 
puted on IBM Relay calculator. 
Argument in grades or units of com- 
plete circle; $4 = 100°; 2x = 1° 
sin x, cos x; 0.25 — x, A: 
= 0(0°0001)0.125; 4D; 5D; and 
6D 
sin x, cos x: x = 0(2001)1; 5D 
First difference given at (1). 
sinx, cosx, A*: x = 0({01)1; 8D 
(third order interpolation) 
sin x, cos x, A, &; 0.25 — x: 
x = 0(£00001)0.125; 15D 
Rounding error in last place may oc- 
casionally be greater than 3} unit. 
At(14) C, A? is given in place of &. 
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Source 
(26) 


[(47)-(18)] 


[(26)-(21) ] 


[—(3)] 
[(47), 
(28)-(28) ] 
[(47), 
(28)-(28) ] 
(—(3)] 
[—(3)] 
[—(3)] 


[(45)-(4) ] 


[(47), 
(28)—(28) ] 
[(80)-(17)] 
[(23)] 
[(47), 
(28)—(28) ] 
[(47), 
(28)-(28) ] 


[(47), 
(28)—(28) ] 
[(47)-(28)] 


[(47)-(70) ] 


[(23)] 
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Available at 
(26) 


(18)C 


(21) 


(24)C 
(28) 


(28) 
(24)C 
(24)C 
(24)C 
(4) 
(28) 


(17) 
(23)C 
(28) 


(28) 


(28) 
(28) 
(23)C, (26), 
(28) 


(23)C 


7.31 


7.35 


7.4 


7.44 





Description of Tables 


sin x, cos x: +x = 0.001(2001)99.999; 
8D 
(First differences for functions at 
intervals of °01) 
sin x, cos x: x = 0(1°)50; 20D 
Trigonometric functions; argument in 
mils 
sin x, cos x, tan x, sec x, csc x, cot x, 
$tanx: x = 0(.1)800 mils; 5D 
(with extension to 1600 mils) 
$ cot x, x = 50.9(.1)800 mils; 5D 
Trigonometric functions; argument in 
degrees, minutes, and seconds 
sin x, cos x, A, &:x = 0°(10’)360°; 5D 
cos x, A per sec.: x = 0(1’)180°; 5D, 
6D 
1 — cosx: x = 90°(1’)180°; 5D 


tan x, cot x; A, 6?:x = 0°(10’)360°; 5S 
sec x, csc x; A, 6: x = 0°(10’)360°; 5S 
sin x, cos x, tan x, cot x, sec x, csc x: 
x = 0°(10’)360°; 5D or 5S 
sin x, cos x, tan x, cot x: 
x = 0°1'(1’)89°59’; 5D or 5S 
tan? x: x = 0(1’)88°; 6D 


tan x: x = 0(10’’)30°; 7D 

cos x, A/60: x = 0°(1’)90°; 8D 

tan x, cot x, A per sec: x = 0(1’)90°; 
8D 

sin nx, nm = 1, 2, 4, 6, 8, A per sec: 
x = 0(1')90°; 12D for m = 1; 10D 
to 7D for n > 2. Also other sets to 
4D, 5D, and 6D for some values 
of n. 

sin? x, cos? x: x = 0(1’)90°; 10D 


sin x, cos x, tan x, cot x, A per sec.: 
x = 0(1')90°; 15D 
sec x, csc x: x = 0°(1’)90°; 10D 


sin y, cosy, A/60; y = (b/a) tan x, 

where a and b are the major and 
minor axes of the ellipse on Clarke’s 
spheroid (1866): x = 0(1’)90°; 8D 








[(i1 


[(23 


[(S6 


[ (24 


[(1 
[(2. 


[i 


[(2 
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Source Available at 


[(11)] (11)C 


[(23)] (23)C 


[(S6)-(21)] (21)C 


[(24)] (24)C 


[(18) J (18)C 
[(23)-(7)]_ (NC, (24)C 


[(18)] (18)C 


[(23)-(18)] (18)C 


7.63 


7.8 


7.9 


8.5 


9.10 


Description of Tables 
sinx cosx d* (==) : 


x ° = *'de 
x = 0(.1)20; 8D 
Area and circumference of circle 
A = 4nD*; C = cD: A, 
D = 0(.0001)0.9999: 8D in A; 7D 
in D 


6 
(a) £(0) = [sect tat 1 — 6 = £/b: 


= 0(1’)87°; 4D to 7D 
(b) sin 0(€), cos 0(—), where 6(£) is the 
inverse of £(6): — = 0(.01)2(.05)50; 
8D or 9D 


(c) X(B, 6) = fi cos {e[b(1 — #)]} a 


(4) (8, b) = f° sin (ofc - 9) F 
(76,0) | 
-f cos {6[b(1 — #) ]}t-tdt 


X, Y, T: B = 0.02(.02)3; 
b = 0.1(.1)2; 8D 











x 








f(cos x, a, k) 
cos x 
= k+cosx Je a : 
2 (1 + a? — 2acos x)! 


cos x = —1(.01)1; 5S or 6S 

a = 0.015, 0.016, 0.017, 0.018, with 
k = 12045.4 

a = 0.0000425, with k = 0.02480289 

Tables with arguments in time. See 
28.5 

Inverse circular functions 

arcsin x: x = 0.005(.01)0.995; 7D 

arcsin x, arcsin (x —h), arcsin (x—2h), 
50”, 6_? 

x = 0(.0001)0.989(.00001)0.99999; 
12D 

x, = arcsin f; 
f = 0(.01)1(—.01) —1(.01)0; 2D 

x2 = arccos f; f = 1(—.01)—1(.01)1; 
2D 


X1y, Xe = 2nd, O<A <1; 2D in A 
(critical table of sin x, cos x) 

2x = arcsin x: x = 0.005(.01)0.995; 
5D ind 
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Source 


[(68)-(14) ] 
[(68)-(14) ] 
[(18)] 


[(23)-(11)] 
[(26) ] 
[(23)-(24)] 


[(23)-(10) ] 


[(23)-(27)] 
[(42)-(24) ] 


[(23)-(24) ] 
[(24)] 
[(21)] 


[(18)] 


[(23)-(22) ] 


[(23)-(13) ] 
[(23)-(13) ] 


[(23)-(13)] 
[(23)-(4)] 


[(45)-(18)] 
[(23)] 
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Available at 
(14)C 


(14)C 


(18)C 9.12 


(11)C 
(26) 
(24)C 


(10) 


(27) 

(23)C, (24)C 

(16), (24)C 

(24)C 

(21) 9.15 


(10), (13), 
(18)C 


(22)C 


9.7 
10 


(13)C 10.0 
(13)C 


(13)C 
(4) 


(18) 
(24)C 
(16), (24) 


Description of Tables 


y = arccosx and arccos —x (y in 
degrees) : x = 0.99899(.00001)1; ac- 
curacy to within 0°.01 

y = arccos x (y in degrees): x = 0 
(opt.int.)0.99899; accuracy to with- 
in 0°.01; 95 cards 

arctan x:x = O(opt.int.)100; 6D; 110 

cards 
arctan given in radians 

arctan x: x = 0(.001)7; 6D 

arctan x: x = 0(.001)1; 10D 

arctan x, V, &: x = 0(.001)7(.01)- 
50(.1)300(1)2000(10)10,000; 12D. 
[V is the “backward” first differ- 
ence. | 

arctan x: x = 0(.001)1(1)10(10)- 
100(100)1000; 12D 

arctan x: x = 0(1)2000(10)9000; 12D 

arctan m/n, arccot m/n, m? + n?: 

m = 1(1)n; nm = 1(1)100; 12D 
arctan x: x = 1(1)2000; 15D 


n arctan 1 = na :m = 0(1)100; 20D 


y = arctan x: x = 0(.001)3.75(.01)- 
18.25; (y in degrees); accuracy to 
0°.01. [Extension for x > 18.25 by 
special table. ] 

y = arctan x:x = O(opt.int.)1000; (y 
in degrees) accuracy 0°.001; 350 
cards 

y = arctan x, A: x = 0(.001)1000; (y 
in degrees); 6D 

General spherical triangle (see 28.4) 

Exponential and hyperbolic functions 

e*, e~*: x = 0(.00001)0.00099; 7D or 
8D 

e*, e-*: x = 0.001(.001)0.999; 7D or 
8D 

e*, e-*: x = 1(1)15; 8S 

e*, e-*: x = 0(.0001)1; 10D for e*; 9D 
for e~* 

e~*, A: x = 0(.01)3(.1)10; 6D 

e~*: x = 0(.1)10; 10D 

e-*: x = 0(.001)2.499; 10D 





So 
[(65) 
[(79) 
[(23) 
[(27) 
[(23) 
[(23) 
[(23) 
[(23) 
[(27) 
(79) 
[(14 
[(26) 
[(23 
[(26 
[(23, 
(18, 
[(23, 


[(24 
[(23, 
[(23, 
[(18 
[(18 


[(18 


[(23 





Source 
[(65)—(24) ] 
[(79)-(4)] 
[(23)-(21)] 
[(27) ] 
[(23)-(14) ] 
[(23)—(24) ] 
[(23)—(14) J 
[(23)-(24) ] 
[(27)-(26) ] 
[(79)—(14) ] 
[(14) ] 
[(26)] 
[(23)-(4) ] 
[(26) 
[(23)-(12)] 
[(18)] 
[(23)-(21) ] 


[(24)] 
[(23)-(4) ] 
[(23)-(21) J 
C(18)] 
[(18)] 


[(18)] 


[(23)-(18) ] 
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Available at 
(16), (24)C 


(4) 
(21) 


(27) 


(14)C 
(24) 
(14)C, 
(21), (24) 
(16), (24) 


(23), (26) 
(14)C 
(14)C, (24) 


(26) 


(4) 
(26) 
(12)C, (24) 


(26) 
(18)C 


(21), (23) 
(16) 

(24) 

(4) 

(21) 
(18)C 


(18)C 


(18)C, (23) 


(16), (18)C 


10.33 


10.4 


11 


11.1 
11.2 


11.26 


Description of Tables 


e*, €-7; A, &, A*: x = 2.5(.001)10; 10D 

e*, €-*: x = 1(1)100; 10S 

e*: x = 2.5(.001)5(.01)10; 15D for 
x <5;12D forx>5 

e*: x = —100(.01)100; 15S. Com- 
puted on IBM Relay Calculator. 

e*, €*: x = 0(10-*)0.0001; 18D 

e-*: x = 0(.01)2.49; 18D 

e*, A, A’, A®, At: x = —2.5(.0001)2.5; 
18D for x < 1;15D forx > 1 

e*: x = 0(.001)2.499; 18D for x < 1; 
15D forx >1 

e*, €-*: x = 0(.01)100; 14S and 18S 

e*, e-*: x = 1(1)100; 19S 

(1 — e-*)/x: 0 < x < 13.0369; 6D; 
(pseudo-opt. int.); 792 cards 

sinh x, cosh x, A, &: 
x = 0(.0001)0.01(.001)1(.01)99.99; 
8D 

sinh x, cosh x: x = 0(.0001)1; 9D 

sinh x, cosh x: x = 0(.01)16.11; 16D 

tanh x, coth x: x = 0(.0001)2; mostly 
8S 


tanh x, A: x = 0(.01)9.9; 8D 

tanh x: x = 0(.005)6; 9D 

Natural logarithms 

log. x: x = 1(1)100,000; 16D 

log. x: x = 2(.01)6.6; 7D 

log. x, A: x = 1(.001)9.999; 7D 

log, x: x = 0(.0001)1; 16D 

log. x: x = 0.0001(.0001)10; 16D 

log, x, (radix table): 
0.001 < x < 9.9999; 5D; 76 cards 

log. (— log. x): x = 0(.0001)0.001- 
(.001)0.999(.0001)0.9999; 4D 

log. (— log, x); A, A?: x = 0(.00001)- 
0.006 (.0001) 0.08 (.001) 0.84(.0001)- 
0.984(.00001)0.99995; 5D 

Inverse hyperbolic functions: 

arctanh x = $ log. [(1 + x)/(1 — x)]; 
arc sinh x = log, (x + Vi + x*) 

arctanh x: x = 0(.002)0.5; 6D 
(A also on cards at 18) 
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Source 


[(S0), 
(14)-(14) ] 


[(S0), 
(14)-(14) ] 
[(10)] 


[(14)] 


[(23), 
(14)-(14) ] 
[(23)-(16) ] 
[(23)-(16) ] 
[(23), 
(2)-(2)] 


[(23), 
(24)—(24) J] 


[(23)] 


[(23)] 
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Available at 
(14)C 


(14)C 


(10) 


(14)C 


(14)C 


(16)C 
(16)C 
(2) 


(24)C 


(23)C 


(23) 


11.8 


13 


13.2 


13.4 
13.5 
13.52 


13.8 


13.9 


Description of Tables 


arcsinh x, arc tanh x: 

x = 0.00001(.00001)0.001; 20D 
and x = 0.001(.0001)0.1(.001)0.999; 
7D to 10D 

arcsinh x, arccosh x 


x = 1(.001)3(.01)10(.1)20(1)50; 7D 
and 8D 
Logarithms of complex numbers 
f = log.2,2 =x+ ty: 
x,y = —4(.02)4; 5D 
Addition logarithms 
log. (1 + e*): x = O(opt.int.)13; 5D; 
270 cards 
Exponential integrals for real and 
complex arguments; sine and cosine 
integrals 
zte 
E,(x) = lim [=f wean; 
ca+0 z 
c real 


E,(x): x = 0.01(.01)10(.1)90; about 
9S 
Si(x), A, #&: x = 0(.01)99.99; 10D 
Ci(x), A, &: x = 0(.01)99.99; 10D 
Modified cosine integral: 
—Ci(x) + 7 + log. x 
x = 0(.001)10(.01)49.99; 6D 
E,(x): x = 0(.001)2; 7D 
x = 2(.001)10; 5S or 6S; 
x = 10(.001)16; 5S to 8S 
Exponential integrals for complex 
arguments; z = x + iy 
E,(z): x = 0(.02)4; 
y = 0(.02)3(.05)10; 6D 


 Ex(z) + log. : x, y = 0(.02)1; 6D 


e*E;(z): x = 3(.1)10; y = 0(.05)10; 
6D 
e*E,(z): —x = 0.5(.5)4.5; 
y = 0(.1)4(.5)10; 6D 
—x = §(.5)10; y = 0(.5)10; 
6D 
x = —10(1)10; 
y = 10(1)20; 6D 
+x = 11(1)20; y = 0(1)20; 
6D Mf. 
[Branch cut on negative imaginary 
axis ] 








Sou 


[(24)_ 


[(63), 
(21)-( 
[(23)_ 


[(23)- 


[(23)- 
[(23)- 


[(23)- 


C(4)] 


[(62) 


[(23) 








Source 


[(24)] 


[(63), 
(21)-(21)] 
[(23)] 


[(23)-(14) J 


[(23)-(14) ] 
[(23)-(14) ] 


[(23)—(14) J 


[(4)] 


[(62)-(14) ] 


[(23)-(14)] 


A GUIDE TO TABLES ON PUNCHED CARDS 203 


Available at 

(24)C 
14 

(21) 14.2 

(23)C 14.6 
15 

(14)C 15.1 

(4), (14)C, 

(24) 

(14)C 

(4), (14)C, 

(23), (24) 

(26) 

(24) 

(4) 15.281 
15.411 

(14)C 


(14)C, (24) 15.521 


Description of Tables 


Available at (23); only partially 
punched at time of publication of 
this guide. 

Ex(z): —x = 0(.1)3.1; y = 0(.1)3.1; 
10D 


E,(z) + log s: —x = 0(.1)1; 
y = 0(.1)1; 10D 
Gamma function 


100 + logie I'(x): x = 1(.01)16.99; 7D 


Log. T(z), 2=x+ty: x = 9(.1)10; 
y = 0(.1)10; 14D or 15D 

[Extension of range to x = 0(.1)10 is 
in progress at time of publication of 
this guide. ] 

Probability functions; Hermite poly- 
nomials; moments 


(x) = ee (—$x*); 
a(x) = f 2(t)dt 


2(x), 4a(x):x = O(opt.int.)4.751; 5D; 
886 cards 

Prepared for linear interpolation; 
fixed interval. 

2(x), 4a(x), A, &: x = 0(.001)7.8; 8D 

At (14)C without differences. 

2(x), a(x): x = 0(.01)7; about 10D 

[with auxiliary functions, prepared 
for third-order interpolation ]. 

2(x), a(x), A: x = 0(.001)7.8; 15D 

At (14)C and (26) without differences. 


ee”: x = 0(.05)6.25; 8D 


A = Be*;B= f exp fdt: 
0 
x = 0.1(.01)9.99; 10D in A; 8S, 9S, 
or 10S in B 
Inverse tables 


x(p), where p = 3[1 + a(x) ], and 2(x) 
p = 0.5(.0001)0.9999; 8D 
ie (e-*): x = 0(.1)5.9; 

m = 1(1)15; (24 — )D 
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Source Available at 
[(23)-(14)] (14)C, (24) 


C(14)] (14)C 


[(14)] (14)C 


[(26)] (26) 


[(14)] (14)C 


C(14)] (14)C, (24) 


[(14), (24)] (14)C 


15.523 


15.6 


15.7 


15.8 


15.9 


17-20 


Neseneaiens of Tables 


Vie ~ (ei); x = 0(.1)8.4; 
n= 5 ee — "aes 


= x(x) and —~ = ste) (Hermite 


a5 5= P 
ag 
n = 1(1)10; x = 0(.01)12; about 6S 
Hermite polynomials: » = 1(1)10, 
x = 0(.01)12; exact 


| ae, 
me) = Gye J, 


X exp (—4f)dt [see (40) ] 
M; = ma ($5 : 5D; exp (M;); 6D; 
exp ag M,); 2D 
[exp (— M;) — exp(— Miz1)]/ 
exp(— M;); 6D 
[M; exp (— Mj) 
— Mizrexp (— Mix:))/ 
exp (— M,) — exp (— Mi41)]; 6D 
1/M = 0.24(.01)0.99; ¢ = 1(1)8 


Fate) =a] f° ao%e-var] pram) 


m = 2(1)10; x = O(opt.int.)44; 5D; 
7637 cards 

Probability integrals relating to com- 
plex arguments 


2(x, y) + ids(x, y) 
= f exp (ity)z(é)dt; 5D 


salz, 9) = f [sin ty e(Oae; SD 
x = 0(.1)4.5, y = 0(.1)0.9; 


x = 0(.05)2(.1)4.5, y = 1(.1)1.9; 


x = 0(.05)3(.1)4.5, y = 2(.1)2.9; 
x = 0(.02)1.5(.05)3(.1)4.5, 


y = 3(.1)3.9; 

x = 0(.02)2(.05)3.5(.1)4.5, 
y = 4(.1)4.9; 

x = 0(.02)2.5(.05)3.5(.1)4.5, 
y = 5(.1)6. 


Altogether 6216 cards. 
q(R, x) — f te +2") T(t) dt 
R 


(Offset circle probabilities) 
R = 0.1(.1)20; x = 0(.05)1; 6D 
Bessel functions 








[(6 


[(6 
[(5) 
[(5) 


[(2. 


[(6 
(16) 


[(6¢ 
[(6¢ 
[(6¢ 
(16) 


[(2: 
[(s: 
[(6¢ 
[(2: 
[(6 
(16) 
[(2: 


[(6¢ 








Source 


[(S)-(16) ] 
[(S)-(16) J 


[(23)] 


[(66)—(14) ] 


[(66)-(11) ] 
[(S)-(16) ] 
[(S)] 


[(23)-(16)] 


[(66), 
(16)—(16) ] 


[(66)-(11)] 
[(66)—(14) ] 
[(66), 
(16)—(16) ] 


[(23)] 
[(S3)-(24) ] 
[(66)—(24) ] 
[(23)-(14) ] 
(16)-(16) ] 
[(23)-(14) ] 


[(66)—-(24) ] 
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Available at 
(16)C 
(16)C 


(14)C, (24) 


(14)C, (24) 


(11)C 
(16)C 
(5)C 


(16)C 


(16) 


(11)C 
(14)C, (24) 
(16) 


(23)C, 
(24)C, (26) 
(24) 

(24)C 


(14)C, (23)C 


(16) 


(14)C 


(24)C 


17.1 


17.2 


17.3 


18.1 


18.2 


18.3 


Description of Tables 


Jo(x), Ji(x), A: x = 0(.1)25; 6D 
Jo(x), Ji(x), ATo(x): x = 0(.1)99.9; 
6D 


Jo(x), Ji(x): x = 0(.01)10; 10D 
First four advancing differences 
given along with functions at 
(14)C. 
Jo(x), Ji(x): x = 10(.01)25; 10D, A, 
A?, A}, A‘ 
Jo(x), Ji(x): x = 0(.01)25; 10D 
J2(x), J3(x): x = 0(.1)25; 8D 
J,(x): 
n = 0(1)15; x = 0(.001)25(.01)99.99; 
18D for n < 3; 10D for > 4 
n = 16(1)100; x = 0(.01)99.99; 10D; 
J,(100): 2 = 0(1)100; 10D 
Note: At the time of preparation of 
this guide, values for m = 79(1)100 
and J,(100) were not yet available 
for distribution. 


AE June): n = 1, 2; x = 0(.1)16; 
8D 

Y,(x), A: m = 0(1)3; x = 0.01(.01)1; 
6D 


[Yo(x) and Yi(x) checked ] 

Yo(x): x = 0(.01)25; 8D 

Yo(x), ¥i(x): x = 10(.01)25; 8D 

Y,(x), A: #2 = 0(1)3; x = 0.5(.1)25; 
8D 

[Yo(x) and Yi(x) checked ] 

Yo(x), Yi(x): x = 0(.01)10; 10D 


Y, (x): x = 0(.1)10; m = 1(1)21; 10S 
Ip(x) and I,(x), #: x = 0(.001)5; 8D 
Io(x), Ii(x): x = 0(.01)10; 10D 
I,,(x) or e~*I,,(x): nm = 0(1)3; 
x = 0(.1)20; 8D 
[e*I,(x) for x > 5; entries for 
n = 0, 1 checked ] 
I(x) and I_4(x); x = 0.01(.01)25; 
10S 


Ko(x), K.i(x)é; also auxiliary func- 
tions for interpolating near the 
origin: x = 0(.01)5; 7D to 10D 
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Source 
[(66), 


(16)—-(16) ] 
[(48)-(21) ] 


[(23), 
(14)-(14) J 


[(14)] 


[(14)] 
[(14) J 
[(66)-(24) ] 


[(66)-(24) ] 
[(24)] 
[(14)] 


[(14)] 


[(23)] 
[(23)] 


-[(23)] 
[(24)) 


[(24)] 
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Available at 
(16) 


(21) 
(14)C 


(24)C 


(24) 
(24) 
(24) 


(24) 
(24) 


(24) 


(24) 


(23)C, 


(24), (26) 
(23)C, (24)C 


(23)C, (24) 


(24) 


(24) 


18.4 


19.0 


20.0 


Description of Tables 


K,(x) or e*K,(x): n = 0(1)3; 

x = 0.5(.1)20; 8D 
[e*K,(x) for x > 5; entries for n = 0 

and 1 checked ] 
e~*Io(x), e~*I;(x), e*Ko, e*Ki, e*: 

x = 0(.02)16; 7D 
e~*Iy(x), e~*I,(x): 

0 <x < 733.4104599; 8D or better; 

1996 cards 
[Quasi-optimum interval, arranged 

for quadratic interpolation, with 

modified first and second divided 
differences. ] 
e-*, e-*I (x), e771, (x); A, &: 

x = 0(.01)2.49; 18D for e~*; 8D for 

other functions 
e~*Io(x): x = 0(.01)36; 8D 
e~*I,(x): x = 0(.01)72; 8D 
e*Io(x), e-*Ii(x), e*Ko(x), e*Ki(x); 

8: x = 5(.01)10(.1)20; 8D 
e~*I (x): x = 2.5(.01)5; 9D 
e~*I,(x): x = 20(.5)120(1)625; 8D to 

10D 
e-*I (x): x = 20(1)48; 37.5(2.5)92.5; 

75(5)175; 160(10)380, 325(25)800; 

10D 
e~*I,(x): x = 16(1)48, 37.5(2.5)92.5; 

70(5)170; 130(10)430, 325(25)700; 

8D, 9D or 10D 
Bessel functions for complex argu- 

ments 
Yo(z) and Yi(z): z = p exp (i¢); 

p = 0(.01)10; @ = 0(5°)90°; 10D 
Yo(ix), Yi(ix): x = 0(.01)10; 10D 
Bessel-Clifford functions 
a J,(2Vx), and x Y,(2Vx): 

n = 0,1; x = 0(.02)1.5(.05)3(.1)- 
13(.2)45(.5)115(1)410; 8D or 9D 
Io(2Vx), [Is(2Vx)]/Vx; #; Ko(2Vx), 

[K,(2Vx)1/Vx: 

x = 0(.02)1.5(.05)6.2 
e-¥2],(2Vx), [e*¥#1,(2Vx) 1/ Vx, 

eM#K (2Vx), [e”#Ki(2Vx) 1/Vx, 3: 
x = 6.2(.1)13(.2)36(.5)115(1)160(5)- 

410; 7D to 9D 








[(49) 


[(78 


[(61 


[(23 


[(72 
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Source Available at 


[(49)-(24)] (24)C 20.8 


22 


[(78)-(16)] (16)C, (24) 22.1 


[(61)-(24)] (24) 22.2 


[(23)] (23)C, (24) 


[(72)] (24) 


Description of Tables 


f Jo(Au) exp tudu 
= Jr, x) + iJ,(r, x) 


f : Yo(Au) exp (iu)du 


= N.(A, x) + iN,(A, x) 
Auxiliary functions: 
CAA, x) = N.(A, x) 


- : [In (Ax) ]J-(A, x) 
C.(A, x) = N,(A, x) 
-- 2 [In (Ax) JJ,(A, x) 


JAR, x), Js(A, x), NCA, x), NCA, x): 
x = 0(.02)2(.1)5; A = 0.1(.1)1; 6D 

Auxiliary functions: 

x = 0(.02)x.; x. < 1.6; 6D 

Riemann Zeta function; 
functions 

(s — 1)f(s);3 = x + ty; x = 0(.05)2; 
y = 0(.05)4; 6D 
[computed on the EDSAC] 

(a) Characteristic values, be,(s) and 
bo,(s), associated with even and odd 
periodic solutions, respectively (pe- 
riod x or 2x), of Mathieu’s equation 

y’ + (6 — scosx)y = 0 
b = be,(s), bo,(s), ®; O << s < 100; 
r < 15; 8D 

(6) Trigonometric coefficients 
De, (s) and Do, (s) associated 
with periodic solutions, same range 
of r and s, but at larger intervals; 
9D or 9S 

(c) Joining factors, relating various 
solutions, 5?; mostly 8S, same range 
of s and r. For precise definitions 
see (61). 

(d) Periodic solutions Se,(s,x) and 
So,(s,x);0 <s < 100,r = 6(1)15; 
x = 0(1°)90°; 7D 

Solutions g(t) and A(t) of 
y”’ + 1 + kcost)y = 0; 

g(0) = h’(0) = 1; h(0) = g’(0) = 0 

e = 1(1)10, = 0.1(.1)1;¢ = 0(.1)3.1; 

3.14(.0004) 3.1428 


Mathieu 
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Source Available at 
[(13)] (13)C, 
(16), (18)C 
[(23)-(4)]  (4)C, (24) 


[(S9)-(15)] (15)C 


[(13)] (13)C 


[(20), (20) 
(23)-(20) ] 


22.73 
23.1 


23.82 


24.2 


Description of Tables 


e = 30(10)90,200, 300, 400; k = 0.02; 
e = 900, k = 0.01 

« = 5, k = 0.12(.02)(.24) and 
1.1(.1)1.9 

Functions tabulated: g(#), g’(t), h(t), 
h’(t); 3D to 5D 

[computed on the ENIAC] 
See 24.2, Harmonic analysis 


4-Pt. Lagrangean interpolation coeffi- 
cients: x = 0(.0001)0.5; 8D 
[7D at (18) ] 
[from symmetry, this is equivalent 
to giving coefficients up tox = 1] 
6-Pt. Lagrangean interpolation coeffi- 
cients: x = 0.001(.001)0.999; 10D 


Orthogonal polynomials for curve 
fitting 
Eo’ (x) = Yo; E:'(x) = Au(x — 2); 
= = 3(n + 1) 
[Orthogonal with respect to sum- 
mation from 1 to n. ] 
£,/(x) and AR: k = 1(1)5;” = 3(1)104; 
exact 
See 23.82 of (40); MTAC, v. 1, p. 
148-150, (59), and (58), p. 307-327, 
for uses of these polynomials. 


Harmonic analysis 
A sin 2rnx, A cos 2rnx 
n = 1(1)30; x = 0(1)120; 
A = 1(1)10(10)100(100)900 
5D for A = 1; 1D for others 
A sin 2xnx, A cos 2rnx 
n = 0(1)30; x = 0(.002)0.25; 
_ 2A = 1(1)5; 10(10)50; 
100(100)500; 3S—5S 
1 \15 
n = 0(1)20;x = (5) : 
+A = 1(1)5; 10(10)50; 100, 200, 
500; 1S-3S 
1 15 


n = 0(1)40; x = 0( 355) a0 


1 
+A = aqnpgs, Ser ane Lee 


1S-3S 
A sin 2xnx: n = 1(1)143; 
x = 0.0005(.0005)0.25; 3S—5S 





C(18 


[(18 


[ (24 


((7! 


C(14 


[(1 


C(14 


[(4) 


[(14 








Source 


[(18)] 


[(18)] 


[(24)] 


[(71)-(12)] 


[(14)] 


[(14)] 


[(14)] 


[(4)] 


[(14)] 
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Available at 
(18)C 


(18)C 


(24)C 


(12)C, (24) 


(4), (14)C, 
(23), (24) 
(24) 


(4), (14)C, 
(24) 


(4) 


(14)C 


25 
25.1 


25.2 


25.3 


25.4 


Description of Tables 


fi = cos 2; fe = sinz: 

z = 0.01 mn +/20; m = 1(1)600; 
n = 1(1)200; 6D 

fs = (sin z)/z; same range of z; 5D. 
Total cards: 40,200 

Note: Same range of z, at double the 
interval in m and n for f; and f, and 
at four times the interval in m and 
double the interval in n fw f2, to 
four decimals, on 7500 cards. 

Additional table of f;, without f; 

and fe, same range of z; 5D; 13,400 
cards 

fi: = A cos (2xmn/k); 
fe = Asin (2xmn/k): 

m=0,1,2;2 =0,1,2;k =3 

m = 0, 1, 2, 3, 4; 2 = 0, 1, 2, 
k=5 

A = 0(.000001)0.000999 and 
—0.499(.001)0.500; 7D 

sin mx, cos mx: n = 1(1)31; 
x = 0(1°)90°; 15D 


S(k, n) = f x* sin (nwx)dx; 


3, 4; 


C(k, n) = f x* cos (nwx)dx 


S(k, n) and C(k, n); k = 0(1)10; 
n = 1(1)100; 10D 

Random numbers 

Digits of 0(1)9 drawn at random, with 
equal probability. 50 random digits 
per card; 20,000 cards 

Random numbers, positive and nega- 
tive; logs of random numbers, 
numbered serially 0(1)6999 

Random Gaussian deviates, x from 
N(0.1). 

Range: —4.417 < x < 4.417; x to 

3D; x* to 6D 

10 deviates per card, 10,000 cards 

Cosines of random angles; 9D 
Cards 0(1)14,566; random angle 6: 

0 <8 < 9.4247 

Correlated Gaussian deviates from 
N(0, 1) 
Xn4i(p) = pta(p) + V(1—p)Es, 








Source Available at 


(6) see (54) 


[(26) ] (26) 


[(66)-(16)] (16)C 
[(16)] (16)C 


[(21)] (21) 


(24) 


[(60)-(19)] (73) 
Ref. 114 


[(60)-(19)] (73) 
Ref. 101 


[(74)-(19)] (19) 


27 
27.1 


27.2 


27.25 


27.3 
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Description of Tables 


where & is a random Gaussian 
deviate, p is the correlation coeffi- 
cient, and xo(p) is an uncorrelated 
Gaussian deviate. p = 0.6, 0.8, 0.9, 
0.95, 0.97, 0.98, 0.99, 0.995; 
n = 0(1)9999; 3D 

Tables relating to aerodynamics 

A wide range of tables on punched 
cards, relating to shock waves and 
other phenomena, for various Mach 
numbers; ‘“‘Busemann”’ coefficients, 
altitude functions, etc., are avail- 
able at (6). See (54). 

Tables relating to number theory 

Powers of primes. 2*: m = 1(1)100; 
exact; [see also 2.3 ] 

Table of primes p: p < 10,000 


Powers of primes for all primes less 
than 1580 p*: k = 1(1)5; exact; 


N 
alsoi1 + ~? + p*; > *; exact, for 
k=0 


N= 2,3,4,5 

Fermat’s quotients and extraordinary 
primes. See description in MTAC, 
v. 5, p. 84-85, [(73), file 117(F)]. 
n, pb, p, «, f, r, R: nm = 2(1)2962 
[computed on the ENIAC] 

The largest prime factor of x: 
x = 1(1)4200 


27.4 Sums of fifth powers of the divisors of n: 


27.5 


m = 1(1)5000; exact 
[See (60)b] 


Ranamujan’s function +r(m); [See 
(60)a] 
(a) r(n); A OIE z {r(m)}?: 


nm = 1(1)2500; N = 10(10)2500; 
exact 
(b) r(p), pa prime: 1000 < p < 2500; 
exact 
(c) |r()|p-"?; 6D 
Stencils for the solution of systems of 
linear congruences modulo 2; 1024 
cards 


. 



































[(74) 


[(18) 
[(S2) 


C(SS) 


[(76) 
(77)- 


[(17) 
[(17) 
[(22) 
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Source Available at Description of Tables 
[(74)-(19)] (73) 27.6 Table of integer solutions of 
Ref. 125 ly —x|<-x 
Range:x < (10°)/9 or y < (10°)/27 
28 Tables relating to astronomy [also 
see 7.44] 
[(18)] (18)C 28.1 Calendar date to Julian date 1853 
(1 mo.) 1950; 92 cards 
[(52)-(18)] (18)C 28.2 [Kepler’s equation } 


E = M —esin E: e = 0(.01)0.38; 
M = 0(1°)360°; accuracy: 0°01; 


14,079 cards 
[(55)-(24)] (24)C 28.3 Tables for rocket and comet orbits 
(Kepler’s equation) 
[See (55) under Sources, for defi- 
nitions ] 
(a) Elliptic orbits: CU), S.(U), 
XAU), A, #: 


U = 10-*n; 2 = 1, 2, 3, 4, and 
n> 5:0.15 < u < 150 
n=0:0.15 < u < 3.15 
(b) Hyperbolic orbits: C,(U), S,(U), 
X,(U), A, &: 
2 =1,2,3,4and2 = 5: 
0.15 < u < 150 
n = 0:0.15 < u < 300 
mostly 7S and 8S 
[(76), (22)C 28.4 General spherical triangles 
(77)-(22)] Arguments: sides a, b; angle C 
Functions tabulated: (90 —c) (degrees 
and min.) and B (degrees) 
Range: a = 0(1°)89°; b = 61°(1°)90° 
to 119° 
C = 0(1°) upper limit 
The upper limits of 6 and C are such 
that c covers the range 0 to 90°. 
Ac along b. 


28.5 Tables with arguments in time 
[1° = 15”; 1™ = 15’) 


[(17) ] (17)C tan x, sec x: x = 0(1*)20™; 7D 
[(17)] (17)C tan x, sec x: x = 0(#1)20™; 8D 
[(22)] (22) 28.6 Navigation Table H. O. 218 


This table tabulates altitudes (—5° to 
+90°) and azimuth (0 to 360°) as 
functions of declination (—30° to 
+30°), local hour angle (0 to 360°), 








212 


Source 
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Available at 


(30) 


[(23)—(28)] (28) 


[(28)] (28) 

(28) 

(28) 
National Bureau of Standards 


Institute for Numerical Analysis 
Los Angeles 24, Calif. 


29 


30 


30.1 


30.2 


30.3 


30.4 


Description of Tables 


and latitude (—90° to +90°). Com- 
plete information about this and 
similar punched card tables is avail- 
able at (22). 


Actuarial tables 

Tables of “(CSO Monetary Values,” 
comprising a wide range of tables 
relating to annuities, etc., at vari- 
ous interest rates, have been de- 
posited with (30); write for more 
complete information. Additional 
tables based upon an interest rate 
of 24% and the use of continu- 
ous functions, prepared by the 
John Hancock Mutual Life In- 
surance Company, can be obtained 
from (30). 

Tables relating to map projections 
and geodesy 

Meridional arc; x = 0(1’)90°; length 
to 0.001 meters 

Tables for the conversion of latitude 
and longitude to grid coordinates 
and vice versa, on the following 
spheroids: 

Hayford, Clarke (1886) and (1880), 

Bessel: all for lat. 0-80 

Everest: lat. 0-45 
All at 1’ intersections 

Tables for polar stereographic grid 
coordinates, north and south polar 
areas, lat. 79°30’ — 90°. 

Tables for the universal transverse 
Mercator projection (for computa- 

‘tion of plane coordinates to 1 cm. 
accuracy) 

Also see 7.44. For other punched 
cards relating to map projections 
write to (28). 

G. BLANCH 
E. C. YOWELL 
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Formulas for Finding the Argument for Which 
a Function Has a Given Derivative 


The problem of finding where a function f(x) that is tabulated at equal 
intervals h has a given derivative (i.e. given f’(x) = f.’ = f'(xe + ph), to 
find p) arises very frequently in numerical computation, especially in the 
location of a maximum or minimum point when f’(xo + ph) = 0. The usual 
methods for finding where the derivative of a function has a certain value 
are either (1) to calculate the derivatives at the tabular points and then 
perform inverse interpolation, (2) to estimate where the derivative has the 
desired value and then subtabulate the first derivative in which connection 
a table by Sauzer! will be found helpful, or (3) to reverse the series in p, 
obtained by differentiating any one of the direct interpolation formulas for 
f(xo + ph). The drawback in (1) is that it involves two separate steps, and 
the drawback in (2) is that it requires much more work than is necessary, 
together with a large amount of continual personal inspection. Thus method 
(3) which gives p in terms of the tabular entries f(x» + mh) seems to be the 
most convenient, and it is recommended by G1BB,? in the form of reversion 
of a series involving differences. This present note gives the coefficients of 
the reversed series for p directly in terms of the tabular entries, so that no 
differences are necessary. However, if differences are already available, one 
might use them conveniently in the formulas given at the end which are 
obtained by reversing the series for the derivative of the Newton-Stirling 
formula using central differences 5)?" and mean-central differences y5,?"*", viz., 
hf’ (xo + ph) = (udo — sudo + soudo® —---) + p50? — Pydot + godo® —---) 

+ PU hubs? — Jude) +--+) + PAB! — debe! +--+) 

+ PA geudd —-++) + p(xbebe —---) +. 
which, after transposition of the constant term in the right member, and 
division by the coefficient of p, assumes the form r = p + sp? + tp? + up' 
+ vp’ + ---. The formulas are for all cases where 3 through 7 points are 
needed in direct interpolation. These formulas are similar to the formulas 
for inverse interpolation by the author;? but these combine the two opera- 
tions of differentiation and inverse interpolation in a simple manner. 

The function f(x) requires m points for direct interpolation when 
the interpolating polynomial which equals f(xo + mh) = f(xm) = fm» 
Xm = Xo + mh, m = — [(m — 1)/2] to [n/2], is of the (nm — 1) degree. 
Corresponding to m = 3(1)7, quantities r, s, t, wu and v are defined below. 
(The user must note carefully that in f,’ = f’(x), the derivative is with 
respect to x, and not p.) 
3-Point 
r = (2hf. + fr— fi)/D,s =t=u=v =0, 
where D = 2(f_1 — 2fo + fi). 

4-Point 
r = (6hf.’ + 2f_1+ 3fo — 6f1 + f)/D, 
$= 3(— fit 3fo — 3fit fr)/Dt=u=v=0, 


where D= 6(f-1 oe 2fo oa fi). 
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5-Point 

r = (24hf.' + 2(— fio + 8f-1 — 8f1 + fe))/D, 
s = 6(— f_e + 2f_-1 — 2fi1 + fe)/D, 
t = 4(f_2 — 4f_1 + 6fo — 4f1 + f2)/D, 
u=v=0, 
where D = 2(— f_2 + 16f_1 — 30fo + 16f1 — fe). 

6-Point 
(120hf.' + 2(— 3f_2 + 30f_1 + 20fo — 60fi + 15f2 — 2fs))/D, 
s = 15(— f_2 — fr + 10fo — 14f1 + 7f2 — fs)/D, 
t = 20(f_2 — 4f_1 + 6fo — 4f1 + fo)/D, 
u = 5(— f_o+ 5f_1 — 10fo + 10f1 — Sfe + fs)/D, 
v=0, 
where D = 10(— f_2 + 16f_1 — 30fo + 16f1 — fe). 


r 


7-Point 
r = (720hf,’ + 12(f_s — 9f_2 + 45f_1 — 45fi + 9f2 — fs))/D, 
s = 45(f_s — 8f_2 + 13f_1 — 13f1 + 8f2 — fs)/D, 
t = 20(— f_s + 12f_2 — 39f_1 + 56fo — 39f1 + 12f2 — fs)/D, 
u = 15(— f_s + 4f_2 — Sf_a+t 5fi — 4fe + fs)/D, 
v = 6(f_s — 6f_2 + 15f_1 — 20fo + 15fi — 6f2 + fs)/D, 
where D = 4(2f_s — 27f-2 + 270f_-1 — 490fo + 270f; — 27f2 + 2f;). 
Then the quantity p is given by 
p=r—rs+ (2s? — t) + (— 5s? + Sst — u) 
+ r5(14s* — 21s% + 32 + 6su — v) 
+ r6(— 4255 + 845% — 28st? — 285°u + 7itu + 7sv) + ---. 


When differences are available it is convenient to define the r, s, ¢, uandv 
according to the scheme below, the formula for p being the same as above. 
However, these formulas using differences, when expressed in terms of the 
tabular entries f,, are not the same as those previously given directly in 
terms of fm. Thus, when differences are available, these formulas using differ- 
ences might serve as a partial check on those formulas without differences, 
or vice versa: 


3-Point 
r = (hf.’ — do) /i?, s =t =u =v=0. 


4-Point 


r = (hfe! — wdo + gud0*)/d0, 
S$ = fub'/ie,t =u =v=0. 
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5-Point 
r = (hf. — ube + Sud0*)/A, 
= $ube/A, 
= 48!/A, u =0 = 0, 
where A = 5? — pydo*. 
6-Point 


r = (hf.’ — wdo + ude? — goudd®)/A, 
= (3udo® — Sudo®)/A, 

t = $50'/A, 

u = grudo®/A, v = 0, 

where A = 6? — pydo'. 

Point 
= (hf. — ubo + gude® — goudo®)/A, 
= (udo? — gudo°)/A, 

t = (350' — gedo*)/A, 
vxube®/A, 
v = r3950°/A, 


where A= 50" sags qxdo! + sedo°. 


u 


H. E. SALzer 

National Bureau of Standards 
Washington 25, D. C. 

This work was sponsored in part by the Office of Air Research. 

1H. E. Sauzer, “Table of coefficients for obtaining the first derivative without differ- 
ences, ” NBS, Applied Mathematics Series No. 2, 1948. 

2D. Grips, A Course in Interpolation and Numerical Integration for the Mathematical 

Laboratory. London, Bell, 1915. 


3H. E. SALZER, “A new formula for inverse interpolation,” Amer. Math. Soc., Buil., v. 
50, 1944, p. 513-516 
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915[F].—K. Fricntt, “‘Statistische Untersuchung iiber die Verteilung von 

Primzahl-Zwillingen,”’ Oster. Akad. Wiss., math.-nat. Kl., Anz., 1950, 

p. 226-232. 

Information on the distribution of twin primes (p, + 2) and quadruplets 
(pb, p + 2, p + 6, p + 8) is given for the first 1020000 natural numbers. The 
information is based on the old table of CHERNAC.' The main table gives the 
number of prime pairs in each of the 1020 chileads 1000” < p < 1000(m + 1) 
for m = 0(1)1019. There is no chilead devoid of prime pairs and only 
three (nm = 688,851,927) with but a single prime pair. The rows of this 
table are added to give the number of prime pairs in each of the 102 myriads 
10000” < » < 10000(m + 1) for m = 0(1)101. These frequencies are, in 
turn, added by tens to give a 10 entry table for each interval of 100000. The 
grand total gives 8168 prime pairs < 10°. 
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As for the quadruplets, there are enumerated in each of the 10 intervals 
of 100000; the total number is 166. On p. 232 the largest quadruplet, p + 8, 
is given for each of the 172 cases below 1020000. 

Previous tables concerning twin primes are mentioned in MTAC, v. 4, 
p. 84. The reviewer has not as yet attempted to reconcile discrepancies 
between these and the tables under review. No doubt some of these are due 
to errata in Chernac. 


 #. £.. 
1L. Caernac, Cribrum Arithmeticum etc., Deventer, 1811. 


916[I, O].—R. E. BEARD, ‘‘Some notes on approximate product integration,” 
Inst. of Actuaries, Jn., v. 73, part II, 1948, 356-403. 


The author studies approximation formulas of the form: 


ff se 0@<x [Eafe ff o@ar +R, ff oar, 


@ < %1, X2, -+*, Xn <b, where the a,, x, are independent of f(x) and 
Se o(x)dx ¥ 0. $(x) is generally a known, tabulated function and f(x) is a 
complicated or empirically given function. One of the more important appli- 
cations of such formulas is the evaluation of continuous single premiums 
in the field of life and disability contingencies where ¢(x) = v* is the dis- 
count factor and f(x) is derived from a mortality or disability table for one 
or more lives. 

The paper is divided into two parts: a theoretical part in which formulas 
are developed for the determination of a,, x, and the remainder term R,, and 
a part containing tables for the a, and x, based on the formulas of the first 
part for low orders (i.e. small 7). 

In the theoretical part the following cases are considered: 


(1) Neither the a, nor the x, are given; R, depends on f®” (x). 
1 
(2) a3 = G2 = +++ =Q, = = R, depends on f+» (x). 


(3) m = 3, a; = a3; R; depends on f(x). 
(4) The x, are given, R, depends on f(x). 


If m; = Sa x'b(x)dx/ Sa’ o(x)dx are the moments of ¢(x) which are 
assumed to exist and to be known, the x, in the cases (1) and (2) are found 
to be the roots of a polynomial of dégree m whose coefficients are determi- 
nants the elements of which are multiples of the moments m,. In the cases 
(1) and (4) the a, are expressed as quotients of two determinants whose 
elements depend in a simple manner on the x, and the m,. Simpson’s rule, 
Weddle’s and Hardy’s formulas are examples of formulas belonging to 
case (4). 

In deriving his formulas for the remainder term R, the author makes 
an application of the mean value theorem which is correct only if $(x) is 
non-negative and 0 < a < 6 and arrives thereby at an expression for R, of 


the form f(é;) = — f™(&)h(x1, «++, x-) where k has the values indicated 


above, h(x1, ---, xr) is a function of the x, only which in many cases can 
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be expressed as a simple function of the m;, and 0 < &;, & < b. He then 
simplifies this expression by stating that, for practical purposes, £, and £2 
may be replaced by a common value with 0 < — < 6 which is evidently 
correct only in particular cases. The remainder terms given by the author 
are therefore useless in many cases for the purpose of measuring the accu- 
racy of the approximation formulas. It is, in any event, doubtful whether 
remainder terms expressed by means of derivatives of high order are useful 
for empirically given functions f(x). 

The author proceeds to develop from the general formulas the formulas 
corresponding to low values of m. In many cases he first makes a linear 
transformation of the variable x to a variable X so as to make m, = 0 and 
m, = 1 and then expresses the remaining moments by means of Pearson’s §;. 
In those formulas where 6; for i > 3 appear he makes the further assump- 
tion that the 6; for i > 3 can be expressed as functions of 6; and 2 in the 
same manner as for Pearsonian distribution functions. 

In another special application the author considers special functions ¢(x) 
such as $(x) = 1, e*, exp (— {(x — m)/c}*), (1 — x)™(1 + x)™. 

No attempt has been made to classify the various formulas by some 
optimum properties or by the magnitude of the remainder terms, the only 
test being a comparison of the true values of some actuarial functions with 
the results of various approximation formulas. 

In the numerical part of the paper the following tables are given: 


Table 1: Solutions to 6D for X1, Xe, a1, a2 for the two point formulas of 
case (1) corresponding to 8; = 0(.1)3.0(.1)3.0. 
Table 2: Solutions to 6D for X1, X2, X3, a1, @2, @3 for the three point for- 
mula of case (1) corresponding to 6; = 0(.2)2 and 6, = 2(.5)6. 
Table 3: Solutions to 6D for X:, X2, X; for the three point formula of 
case (2) corresponding to 6; = 0(.01).50. 


1 
Table 4: Solutions to 6D for X1, X2, Xs, a1 = a3 = 244 


for case (3) corresponding to 8; = 0(.2)2 and #, = 2.0(.5)6.0. 
Table 5: Moments m, m2, m3, ¢ to 3D and VB: to 5D for the continuous 
function (1 + 4)-* for ¢ = .02(.01).06 and m = 5(5)60. 
Table 6: A similar table for discrete moments. 
Table 7: Solutions to 5D for a1, a2, a3 for the three point formula of case 


(4), where x; = 0, x2 = > » Xm = mand ¢(x) = (1 + 4) for m = 5(5)50 


and «4 = .02(.01).06. 

Table 8: The values to 3D of @, Gz, @szz, Exezz, Szzzzz for the A 1924-29 
ultimate mortality tables for x = 15(1)80. 

Table 9: « to 3D, (8:)? to 5D, 8: to 4D for the function (yl)... for 1 to 
5 lives of the A 1924-29 ultimate mortality table. 


4 
2+A 


» 2 = 


There are also a few auxiliary tables for the purpose of evaluating re- 
mainder terms. 


STEFAN PETERS 
University of California 
Berkeley, California 
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917[I]._Leroy F. Meyers & ARTHUR SRD, ‘“‘Best approximate integration 

formulas,”’ Jn. Math. Phys., v. 29, 1950, p. 118-123. 

For each A = cox(0) + ciux(1) +--+ + Cnx(m) that is an approxima- 
tion to fo"x(t) dt which is exact whenever x(t) is a polynomial in ¢ of degree n, 
m > 1, > 0, there is a kernel function R(¢) such that, when x(é) is of 
class C*+!, 


Rie} = [" xoa— a = [ xemeewar 
0 
The kernel function is defined explicitly by 
R(t’) = ROW] = — Rio] 


where 


a ra 0 if _ft—t)/n! ift<t 
ve = wet) = ty — t')"/n!, ov = dr(t) = { 0 ¥ iff > #. 


By Schwarz’s inequality 


|REx]| < ( f . ic at) ( f c xioto(prat)). 


The best approximation A, for given m and 1, is defined as that which 
minimizes J = fo™k?(t) dt. 

In the present paper, the authors give the best integration formulas for 
n= 1, m = 1(1)20; » = 2, m = 2(1)12; n = 3, m = 2(1)9. For these 
values of m and m, they tabulate the exact values of ¢, ¢1, «++, Cm and J. 
In an earlier paper by Sarp! the best integration formulas are given for 
n = 0, all m, and m = 1, 2,3,m < 6. 

The authors derive some fundamental algebraic relationships between 
the c,’s and J. Then for the case m = 1, recursive relations are derived 
which afford a complete characterization of the best integration formulas 
for any m. Finally, the authors give some conjectures about the con- 
vergence of the coefficients, some of which are true for m = 0 and 1, but 
which are open questions for m > 2. 

H. E. SALzER 
NBSCL 


1A. Sarp, “Best approximate integral formulas; best approximation formulas,”’ A mer. 
Jn. Math., v. 71, 1949, p. 80-91. 


918{I].—Leroy F. Meyers & ARTHUR SARD, ‘‘Best interpolation formulas,”’ 
Jn. Math. Phys., v. 29, 1950, p. 198-206. 


The tabular values x(0), x(1), x(2) of a function x(t) being given, the 
problem discussed here is the approximation of x(u) by an expression of the 
form A = aox(0) + ayx(1) + --+ + aGnx(m), where the value of m and the 
coefficients a = do(u), a1 = a;(u), +++, Gm = Gm(u) are to be determined. 
If A is an exact approximation whenever x(t) is a polynomial of degree n, 
there is a kernel function k(t, wu), such that when x(¢) is of class C**, 


R[x] = x(u) -—A = f xt) (t)k(t, u) dt, 
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where K is the smallest interval containing u and those values of 0, 1, ---, m 
for which the corresponding dp, a1, ---, @m are not zero. For each u, k(t, u) 
is a broken polynomial in ¢ consisting of at most m + 3 arcs, which is defined 


explicitly by 
k(t’, u) = Rov ] = — REoe], 
where 
ie 0 mH ts (¢—#’)*/n! ift<, 
We =e) = a — t')"/n!, ov = de(t) = 0 ift>?t’. 


By Schwarz’s inequality 


4 
IRC]! < a fx corar/ix\ | 


where the modulus M is defined by M = M(u) = [|K| Sxk(t, u)*dt]*. For 
given m, n, u, that A is called best which minimizes the modulus M, and 
it is denoted by Am,n,u- The authors report that for 2 = 0, all m, and n = 1, 
m = 1(1)4, m = 2, m = 2, the conventional polynomial interpolation is 
best, but not for 2 = 2, m = 3,4, 5. 

The determination of Am,x,. involves the prior determination of a num- 
ber of approximations B,,,,,. (for different values of m and u) where Bp, », « 
is that A which minimizes J = fmz{“> k(t, u)*dt. The modulus of B,, .,. is 
denoted by Mm, a, u- 

The authors tabulate the auxiliary function u» = Mz,2./0(u) where 
O(u) = (u — [u])2(1 — u + [wv ])?/120 for m = 2(1)5, uw = 0(.1)2;m = 3(1)5, 
u = 2.1(.1)2.5;2D. The table of » enables one (a) to compare moduli, 
(b) to identify Am+;,2. with B, 2.4, 0 <5 -- m, where each different value 
of m corresponds to a different range of u, and (c) to find the proper argu- 
ment by translation in using the main table. 

The principal table in the article is the collection of formulas for B,,, 2 ., 
0<u<[(m+1)/2], m=2,3,4, and Miou, u < [(m + 1)/2], 
m = 2,3, 4,5. In the Bu2., the coefficients of xo, x1, -+-,%m are given as 
exact polynomials in u. The M2... are expressed as either exact poly- 
nomials in u, or as exact polynomials in u multiplied by @(u). The expres- 
sions for Bn, and M2, >, are different for « lying within different ranges. 

The rest of the paper is concerned with the derivation of those formulas 
and their transformation under a linear transformation of the ¢-axis, 
* = bt+c,b +0. 

H. E. SALzER 
NBSCL 


919[J].—T. M. Cuerry. “Summation of slowly convergent series,’’ Camb. 
Phil. Soc., Proc., v. 46, 1950, p. 436-449. 


This paper studies two transformations of the remainder that are helpful 
in numerical summations. 


The power series considered is >> C,t’, and, in the remainder >> C,t’, 


0 n 
C, is written as a product ¢,f(r), where it is supposed that f(z) has an 
asymptotic expansion 


f(z) ~ Ast + Ays™ + Ags + --- 


ra 


oo ______________ 
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with appropriate conditions. Then, if 


Cn + Capit + Cniol? + --- = ¢n(t) 
and 
d d d 
v=tF Di =F Di=7F) Af(n) = f(n + 1) — f(n), 


the transformations are 
p—1 


: cast fln +1) =X Pbalt)-Da'fln)/1! + Re» 


x Carl f(n +7) = i” = rDe n(t)-A*f(n)/r! + Sa,» 


where Ry, p, Sa,» are ‘‘error-terms.” 

The assumptior of an asymptotic expansion for f(z) is needed for the 
theoretical development, although its coefficients do not occur in the 
numerical applications. On the other hand, although ¢,(¢) is not, in theory, 
very restricted, it must be easily accessible numerically which means, 
virtually, that it must be chosen in such a way that a closed expression 
for it is known. 

A thorough discussion of the error terms R,, », S,,» is given, and upper 
bounds for them are obtained. In practice these upper bounds exceed the 
true error considerably, so that a closer rough estimate has been sought 
which is more useful practically, although it gives the order of the error 
only; for the numerical examples tested the true error rarely exceeds the 
estimate, and never by a factor exceeding about 1.2. 

Two particular families of functions ¢,(¢) are studied in considerable 
detail, and are used in the numerical examples. These are binomial re- 
mainder functions. 

The first family is defined by 
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Bal) = Sti) tG@timtat 77 = (lh — )Bunl?) 
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B,,,(t) = n(n + 1)(n + 2) + (m + 1)(m + 2)(n + 3) oe 
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and"so on. 
The second family, with 





2s—-1 tin-1 


1 3 
bn = 24 ro Qn 2n bam 








an 


for 
exa 


anc 


Uni 


920 








RECENT MATHEMATICAL TABLES 221 











is given by 
B, »(t) = bn + Dasat + Dn+2f? + sited 
=f"{(i-#7-1-—}-—--- —b,77"} 
bn Dasit _2 ‘agns 
B; .(t) = stitsaat 7 2 Sih = @ t)B,.(t)} 
= bn Dasit 
Bu) = GE DmtD | @t+2m+3) 1 
7 b, 
and so on. 


Tables are given of B; »(#) = biset*? = R(r, 0) + il(r, 0) to 4 decimals 


for r = 0.7(0.05)1, @ = 0°(5°)90°, and have been used for the numerical 
examples, which concern the Kapteyn series 


E 2J.(ry) 


r=1 
and are very fully considered. 


J. C. P. MILLER 
University Mathematical Laboratory 
Cambridge, England. 


920[K, L].—B. V. GNEDENKO, Kurs Teorii Veroiatnostei [A Course in the 
Theory of Probabilities]. Moscow and Leningrad, 1950, 388 p. 
On p. 372-385 there are tables of 


o(x) = = #12, x = [0(.01)3.99; 4D] 


(x) = = f , e*2dz, x = 0(.01)2(.02)3(:2)4(.5)5; 4D up to 2.98, 
5D to 8D thereafter 


P,(a) = a a = A(.1)1(1)9, & = [0(1)27; 6D] 


y we a = .1(.1)1(1)3, & = [0(1)15; 6D] 

a0 ! 

P(x) = 5, f" sled, x = 1(1)30, k = [1(1)29; 4D] 
zen ( ) - 





S(x) = (3 ) 1+ ds 
pL +555) 


((n — ayer (2+ 
= 2(1)20, ©, x = [0(.1)6; 3D], x = ©;5D 
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K(x) = © (—1)%e-™, x = .28(,01)2.50(.05)3; mostly 6D 
k=—a@ 





a = .001, .01(.01).05, .1, .15 


6B = [.001, .01(.01).05, .1, .15; 3D] 
R. C. ARCHIBALD 


In 


1-8 
a 


Brown University 
Providence, Rhode Island 


921[L].—AtrcRAFT RADIATION SysTEMS LABORATORY, “Tables of modified 
cosine-integral,”” Stanford Research Institute, 1951, viii + 56 p. 
These tables contain 6D values of 


Ci (x) = fas —28t » 
0 t 


for x = 0(.001)10(.01)50. The computation was carried out, largely on 
IBM machines, by the Telecomputing Corporation of Burbank, California. 
The introduction, by C. T. Tar, discusses the connection of Ci with the 
cosine-integral function and the application of the tables, and describes the 
computation and preparation of the tables. A bibliography is appended. 
The work was sponsored by the U. S. Air Force. 
A. E. 


922[L, M, Q, S].—M. P. Barnett & C. A. Coutson, “The evaluation of 
integrals occurring in the theory of molecular structure. Parts I and 
II,” Roy. Soc. London, Philos. Trans., v. 243A, 1951, p. 221-249. 
Table 1, p. 233. Modified Bessel functions of the first kind, I,4;(x), 

to 7S for m = — 1(1)4 and x = .5(.5)10. 

Table 2, p. 233. Modified Bessel functions of the third kind, K,4;(x), 
to 7S for nm = — 1(1)4 and x = .5(.5)10(1)25. 

These tables are used for the numerical computation of the functions ¢ 
defined by the expansion 


gle? = (tr)-* > (2n + 1)P,,(cos 0)Fm, n(1, t; 7) 


where r? = # + 7? — 2tr cos 6. With these functions the authors form 


Zant t) = [ema ts reat 


and discuss the computation of Z by both numerical integration and ana- 
lytical methods. 

In the memoir it is shown that a large number of integrals occurring 
both in nuclear physics and astrophysics can be reduced to known integrals 
and to Z integrals. Formulas are listed for more than 180 integrals. 

A. E. 


923[L].—C. L. BARTBERGER, ‘“The magnetic field of a plane circular loop,” 
Jn. Appl. Phys., v. 21, 1950, p. 1108-1114. 














RECENT MATHEMATICAL TABLES 223 


The integrals 
i= =f (1 — bcos 6)-'d0 
0 


I; = > f (1 — bcos 6)—! cos 6d@ 
0 


are expressed in terms of complete elliptic integrals. Series expansions are 
also given in ascending powers of 6 and 1 — b. 


Table I (p. 1110-1111) gives 6D values of J;, and table II (p. 1111-1112) 
6D values of I2, for b = 0(.001).809. 
Table III (p. 1113-1114) gives 6D values of J;, I2, I, — In, (1 — bhi, 
(1 — 6b)I; for b = .8(.001)1, and table IV (p. 1114) 6D values of the same 
functions as table III for b = .995(.0001)1. 
A. E. 


924{L].—A. FLETCHER, “Tables of two integrals and of Spielrein’s inductance 
function,” Quart. Jn. Mech. Appl. Math., v. 4, 1951, p. 223-235. 


The tables (p. 229-232) are of 
1 1 
I= f (K — E)dk, J= f (K — E)k*dk 


and —16”(I — a®J)/[3(1 — a)*]. 
Values are given to 10D, 10D and 6D respectively. The range of a is 
0(.01)1. Half a dozen small auxiliary tables are also given. 


925[L].—Harvard University, CompuTATION LABORATORY, Annals, v. 14: 
Tables of the Bessel Functions of the First Kind of Orders Seventy-Nine 
through One Hundred Thirty-Five. Cambridge, Mass., Harvard Univer- 
sity Press, 1951, viii, 614 p. 19.5 X 26.7 cm. $8.00. 

This is the twelfth and final volume of the monumental set of Tables of 
Bessel Functions of the First Order, published by Harvard during the past 
five years—six volumes in 1947, three in 1948, two in 1949 and one in 1951. 
The tabular parts of the volumes fill 7652 pages. The previous 11 volumes 
have been reviewed in MTAC-: v. 2, p. 261-262, 344; v. 3, p. 102, 185-186, 
367, 474-475; v. 4, p. 22, 92. Roughly speaking we now have here 10D 
tables of all J,(x), for x = 0(.01)100, when m = 0(1)111; and for x = 0(.1)- 
100, when has any positive integral value > 111; for m > 135 the values 
of J,(x) are always less than 10-'°. In addition to what is thus stated, for 
n = 0(1)3, x = [0(.001)25(.01)100; 18D]; for m = 4(1)15, x = [0(.001)- 
25(.01)100; 10D]. Detailed information concerning interpolation in the 
whole range is given in Annals, v. 3 and 5; for 10D interpolation the work 
is not excessive. 

Further, in the present volume we have 10D tables of J,(m) for m = 0(1)- 
100. Zero values are not given, since the values found by the Computation 
Laboratory were presented to the Royal Society Committee for use in 
connection with their second volume of Bessel function tables. For »>92, 
x <100, J,(x) +0. For J,(), in the Harvard range, we had earlier: MEISSEL 
(1891), 2 = 20 to 20D; MEIssEt (1895), m = 1(1)24 to 18D; Arrey (1916), 
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nm = 1(1)50(5)100 to 6D; Watson (1922), m = 1(1)50 to 7D; and Havasu 
(1930), = 2 to 101D, m = 10 to 61D, = 20 to 41D, n = 30 to 35D, 
n = 40 to 35D, m = 50 to 30D, m = 100 to 18D. Hence most of the values 
in this special Harvard table are new. 

In the recent Russian table of FADDEEVA and GAvuRIN, RMT 852, the 
argument extends to 124.9, at interval .1, so that some 6D of J,(x) for all 
orders, m = 0(1)120 supplement values given in the Harvard tables. So also 
for 5D zeros < 125, of J,(x); the last zero is for Jis(x). 

The remarkable Automatic Sequence-Controlled Calculator on which 
these tables were computed carried the values of J,(x), for = 0(1)3 to 
23D, and for 2 > 3 to not less than 13D and most of the time much more 
than this; its electromagnetic typewriters also wrote out, for checking 
purposes, 10 differences in every case, and also finally produced the 18D 
or 10D copy which could be sent directly to the printer for offset repro- 
duction. The computation of these tables was only a tiny fraction of the 
work achieved in the ASCC since its activities began in 1945, and have 
continued to the present, 24 hours a day, 7 days a week. 

Only one tabular slip has ever been found in the published volumes, 
J;(72.10) [MTAC, v. 3, p. 41], but this slip was due to some reproduction 
difficulty, and not to any error in computation or in automatic mechanical 
checking. In the last line of the second page of the “Preface,”’ of the volume 
under review, for J(x), m = 0(1)120, x = [0(.01)14.99; 8D], read J,(x), 
n = 0(1)120, x = [0(.1)124.9; 6D]; m = 0(1)13, x = [0(.01)14.99; 8D]. 

In the first page of the Preface, line — 2, for Claire, read Clare. 

These Harvard Bessel Function tables, with most of the values new, 
constitute an outstanding contribution to scientific research. 


R. C. ARCHIBALD 


Brown University 
Providence, Rhode Island 


926{L, S]—Hetmar Krupp, “Bestimmung der allgemeinen Lésung der 
Schrédinger-Gleichung fiir Coulomb-Potential,’’ Akad. Wiss. Leipzig, 
mat.-phys. Kl., Berichte, v. 97, 1950, no. 8, p. 1-28. 


The Schrédinger equation for the Coulomb potential is 
@R _ 2dR |- a3 1 _ WFD |R=0, | 


teat? 2n? 2r? 


dr r dr = @,1,2,-- 


Two solutions are written in the form 


I 
i h0)-= ae aed, $= 4d 


where a = 1+1—2n, 6 = 214+ 2, x = 2r/n. 
With the abbreviations 


(a)o = 1, (a)m = a(a + 1) --+ (@ + m — 1) for m = 1,2, --- 
A = ¥(—a) ifa = 0, —1, —2,--- 
A=V(a-1)ifa#b—1,b-—2,-:- 
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the definitions of the M are 





M(a,b,2) = © Oe= 


x 2M(a, b,x) = [Inx + C — ¥(b — 1) + A]1M(a, }, x) 


= (a).x* = i = tity & 
+E R= =( b+k $3) 











b—2 —)» 
+> , Tia —b+m-+ 1)b'(b — m — 2)ix™-, 
m=o m!T'(a) 


Tables are given for :R, di:R/dr, 2»R, d2R/dr for 0 < x < 15 (the interval 


. . 1/1\7 3/1\7 
in most cases is 5), for 1= 0 and n= 3(4)7,1=1 and n= 3(3)2, 
7 


1 
1=2 anda = 3(3)3 . Graphs of these functions are added. A misprint 


occurs on page 14, line 3 from the bottom and on page 15, line 6; in both 
~etaetit My should be replaced by (" Mi ‘) . However 


K 


the formulas to which the author refers as the source of computation do not 
contain this misprint. 


places ( 


Maria WEBER 
California Institute of Technology 
Pasadena, California 


927[L].—U. E. Kruse & N. F. Ramsey, “The integral fo” y’ exp (— ° 
+ ix/y)dy,” Jn. Math. Phys., v. 30, 1951, p. 40-43. 


In several problems of theoretical physics there arise integrals which can 
be reduced to the real part, I(x), or the imaginary part, K(x), of the integral 
mentioned in the title. Table 2 of this paper gives 5D values of both these 
functions for x = 0(.1).4(.2)8(.5)20. 

Convergent expansions, useful for x < 3, have been given by ZAuN,} 
and Table 1 of the present paper gives numerical values, to 7—9S, of the 
coefficients up to that of x‘. For larger x it is more convenient to use asymp- 
totic series developed by LAporTE,? and 6D values of the first six coefficients 
in these series are also given in the present paper. [The author remarks 
that in TorRREY’s paper* the recurrence formula for the coefficients contains a 
misprint, but Torrey’s numerical values are in agreement with the author's. ] 

Table 2 was computed from these expansions, and V. E. CULLER assisted 
in the computation. 

A. E. 

aS T. Zann, “Absorption coefficients for thermal neutrons,” Phys. Rev., v. 52, 1937, 
. “ 0. Laronre, “Absorption coefficients for thermal neutrons,” Phys. Rev., v. 52, 1937, 
p. — 1%, 


3H. C. Torrey, “Notes on intensities of radio frequency spectra,”’ Phys. Rev., v. 59, 
1941, p. 293-299. 
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928[L, S].—N. MEtTrROPOLIs & J. R. REt7z, ‘‘Solutions of the Fermi-Thomas- 
Dirac equation,”’ Jn. Chem. Phys., v. 19, 1951, p. 555-573. 


Solutions y are given of the equation 


= x(e + pix-t)® (2 = 3-2-4.9-2Z-*) 
in terms of the variable w = (2x)4. The tables are at intervals of .08 and 
extend until y becomes negative. The parameter Z takes on the 24 values 


Z = 6 (4) 14, 16, 18 (4) 26, 29 (4) 81, 84 (4) 92 


and there are 8 different initial slopes y’. Actually 2y is tabulated to 5D. 
The calculations were done on the ENIAC. 
D. H. L. 


929[L].—L. PRANDTL, with the assistance of F. VANDREY, “Flieszgesetze 
normalzaher Stoffe im Rohr,” Zeit. angew. Math. Mech., v. 30, 1950, 
p. 169-174. 
Table 1. 4S table of 


g(a) = cosh a — 5 sinh a + - (cosh a — 1) 


for a = 0(.1)5(.2)10. 
Table 2. 3D table of y(at)/ (a) for = 0(.2).6(.1)1 anda = 1(1)10(2)14. 
Some of the values were obtained by interpolation: these are put in paren- 


theses. 
A. E, 


930[L].—S. Sitver & W. K. SaunpeErs, “The radiation from a transverse 
rectangular slot in a circular cylinder,” Jn. Appl. Phys., v. 21, 1950, 
745-749. 

Table I gives values of ka, n, and ¢o for which Hy’ (ka)/H,’ (ka) < .0001 

and (n¢o)— sin nd > .9. 

Table II gives 4D values of 


° €n0" COS NO €nt” 
ao sin 0H,’ (ka sin 0)/ s=o H,’ (ka) 








for ka = .8, 0 = 10°(10°)90°, 6 = 0°(10°)180°. Here & = 1,€, = 2ifn > 0. 
Table III is similar to table II except that ka = 2.5. 
A. E. 


931[(L].—E. Wo r, “Light distribution near focus in an error-free diffraction 
image,” Roy. Soc. London, Proc., v. 204A, 1951, 533-548. 


In the course of the work the function 


Qan(0) = E (- 1)*LJ.(0) Jom-a(0) + Sega(0) Jamar-e(0) ] 


is introduced, where the J are Bessel functions of the first kind. 
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Table 1 (p. 541) gives 5D values of Q2,(v) for » = 0(1)15 and m = 0(1) 
M(v) where 


M(v) = v +1 forv < 5, M(6) = 6, M(7) = M(8) = 7, M(9) = 8, 
M(10) = M(11) = 9, M(12) = 10, M(13) = M(14) = 11, 
and M(15) = 12. 


The results of some numerical computations involving these functions 
are given in form of diagrams. 


A. E. 


932[V].—E. Gruscuwitz, Calcul approché de la couche limite laminaire en 
écoulement compressible sur une paroi conductrice de la chaleur. Office 
National d’Etudes et de Recherches Aeronautiques (O.N.E.R.A.). Pub- 
lication No. 47, Chatillon-sous-Bagneux, Seine. 1950, 39 p. 

As a generalization of the KARMAN-PoHLHAUSEN' method to steady 
compressible flow, the velocity component wu parallel to a wall y = 0 is ap- 
proximated within the boundary layer by a quartic in » = const fo” pdy 
for a fixed x; moreover, the density p is expressed as a rational function of 7. 
These assumptions lead to the system: 


(1) (0u./v.)d0/dx = F\(K) — (K/bo)[2 — M?F2(k)], 
(2) bo = (1 + .2025 M?)[1 + M?F(K)1/[1 + M?Fi(R)], 
(3) K= (#Pbo/v.) du,/dx, 


with known initial conditions at u, = 0 for the determination of K (which 
is a quintic in Pohlhausen’s parameter A), by) = [p./p_]_~o, and the momentum 
loss @ = JSo"*(pu/p.ue)(1 — u/ue)dy. The velocity u,., kinematic viscosity »,, 
Mach number M, and density p, are supposed known in the main stream 
as well as y, = [y ]um.o0u,, While the functions F;(K), which depend on the 
Prandtl number Pr, are tabulated to 3D (4 = 2,3,4) or 4D (¢ = 1) for 
K = .094(—.001) — .156 [Pr = .725 (¢ = 1, 2, 3, 4) and Pr = 1 (4 = 2)]. 
In addition the displacement 6* = fo" (1 — pu/p.u.)dy of the main stream 
from the wall is given by 6* = 6[bofi(K) + M2f2(K)], the f(K) being 
tabulated to 3D for the above range in K [Pr = .725]. 

In the appendix E. A. EICHELBRENNER describes an exact method and 
presents graphical results indicating that certain aerodynamical quantities 
(excluding the temperature) are equally well approximated by using 
Pr = .725 or the simpler value Pr = 1 [F3(K) = F,(K)] in Gruschwitz’s 
method. 

The tables on pages 13-16 have columns of values of K (as above), 
d (generally to 4S), Fs, Fs, fi, fe, Fi, Fo, Fe [Pr = 1]. There are also rows 
of values for A = + 12, 7.0523, K = — .157 (A, Fi, Fo, F2 [Pr = 1]). 

The identities F, = .595 + F; and fz = .405 — F; should hold through- 
out the table. 

R. R. REYNOLDs 
National Bureau of Standards 
Institute for Numerical Analysis 
Los Angeles 24, California 


1K. PoHLHAUSEN, ‘Zur naherungsweisen Integration der Differentialgleichungen der 
laminaren Grenzschicht.” Zeit. angew. Math. Mech., v. 1, 1921, p. 252-268. 
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193.—CrépiItT COMMUNAL DE BELGIQUE, Tables d’Intéréts et d’ Annuités. 
Brussels, 1950. 
Table V, p. 110, 1/a@j@ .08 for 0.08328680 read 0.08328684. 
L. R. PACKER 
15 Brookland Hill 
London N.W. 11 


England 


194.—J. W. GLoverR, Tables of Applied Mathematics in Finance, Insurance 
and Statistics. Ann Arbor, 1930. 
Part III, p. 394-413, Table of mth derivative of (27)-* exp (— 4). 


p. t n for read 
399 1.14 8 — 30.15397 — 30.154014 
403 2.16 8 13.42550 13.425483 
403 2.44 8 8.04225 8.042232 
407 3.32 7 0.85568 0.855664 
409 3.70 8 — 1.27557 — 1.275600 
413 4.60 8 0.20725 0.207281 
413 4.80 8 0.16334 0.163369 
W. L. SEMon 


Computation Laboratory 
Harvard University 
Cambridge 38, Mass. 


195.—N. R. JORGENSEN, Unterségelser over Frekvensflader og Korrelation. 
Copenhagen, 1916. 
A recomputation of a part of J¢rgensen’s tables of the mth derivative of 
(2x)-* exp (— }x*) by the Automatic Sequence Controlled Calculator re- 
veals the following errata. 


x n for read 
04 6 —5.9506639 —5.95066324 
10 6 —5.7762539 —5.77625459 
14 6 —5.5796134 —5.57961394 
16 6 —5.4581549 —5.45815434 
32 6 — 3.9980943 —3.99809459 
34 6 —3.7642068 — 3.76420646 
40 6 —3.0122141 —3.01241439 
80 6 2.2938198 2.29381943 
82 5 —2.0400226 —2.04002227 
88 6 3.0653608 3.06536044 
1.00 6 3.8715312 3.87153159 
1.02 5 — 1.3734687 — 1.37346845 
1.02 6 3.9619256 3.96192477 
1.06 6 4.1043170 4.10431753 
1.10 6 4.1958467 4.19584621 
1.20 6 4.2103399 4.21033893 
1.22 6 4.1785340 4.17853304 
1.26 6 4.0829542 4.08295339 
1.30 6 3.9475287 3.94752847 
1.32 6 3.8660082 3.86600921 
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for 
3.7759349 
3.4595341 
3.3404608 
3.0852220 
2.6680571 
0.5802503 
2.0712481 
1.0095367 
0.9520677 
0.9679026 
—0.4691537 
0.5801438 
0.9911308 
1.0059216 
—0.0760477 
—0.3491837 
1.0033660 
0.9986367 
0.9917359 
0.9827682 
0.9718380 
—0.2128736 
0.9283347 
—0.8480016 
—0.9216997 
—0.9898676 
— 1.0525184 
—0.1365916 
0.8494399 
— 1.1096856 
0.8267198 
—1.1614158 
—1.2077754 
— 1.2847396 
0.7270882 
—1.3155509 
0.7005090 
— 1.3414083 
— 1.3624566 
—0.1491983 
—0.0214149 
0.6460417 
— 1.3788340 
0.5904438 
— 1.3982379 
0.5624388 
— 1.4015997 
0.5342418 
— 1.4009731 
0.5064238 
— 1.3965440 
— 1.3885091 
0.4509079 
— 1.3770629 
0.4235077 
— 1.3623946 
— 1.3447081 
—1.2755618 
0.2924724 


read 


3.77593384 
3.45953334 
3.34046151 
3.08522283 
2.66805791 
0.58025050 
2.07124870 
1.00953632 
0.95206724 
0.96790228 
—0.46915342 
0.58014344 
0.99113044 
1.00592110 
—0.07604872 
—0.34918347 
1.00336536 
0.99863612 
0.99173665 
0.98276891 
0.97183739 
—0.21287344 
0.92833416 
—0.84800114 
—0.92169927 
—0.98986749 
— 1.05251861 
—0.13659143 
0.84943889 
— 1.10968436 
0.82671889 
—1.16141445 
— 1.20777570 
— 1.28473822 
0.72708742 
— 1.31554946 
0.70050969 
— 1.34140970 
— 1.36245589 
—0.14919850 
—0.02141240 
0.64604256 
— 1.37883586 
0.59044322 
— 1.39823660 
0.56243808 
— 1.40159796 
0.53440589 
— 1.40097219 
0.50642453 
— 1.39654584 
— 1.38851010 
0.45090689 
—1.37705991 
0.42350716 
— 1.36239298 
— 1.34470892 
— 1.27556009 
0.29247276 
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x n for read 
2.56 6 —1.2478101 — 1.24781146 
2.58 6 —1.2180421 — 1.21804283 
2.60 5 0.2437638 0.24376322 
2.60 6 —1.1864509 — 1.18644823 
2.62 6 —1.1532160 — 1.15321833 
2.64 6 —1.1185374 —1.11853985 
2.66 5 0.1756303 0.17563083 
2.66 6 —1.0825918 — 1.08259508 
2.68 3 —0.1232624 —0.12326281 
2.68 5 0.1543471 0.15434760 
2.68 6 —1.0455579 — 1.04556138 
2.70 3 —0.1207053 —0.12070568 
2.70 6 — 1.0076074 —1.00761072 
2.72 5 0.1140524 0.11404817 
2.84 6 —0.7294393 —0.72943953 
W. L. SEMon 
Computation Laboratory 
Harvard University 


Cambridge 38, Mass. 


196.—J. Peters, Eight-Place Table of Trigonomeiric Functions for every 
Sexagesimal Second of the Quadrant. Berlin, 1939; Ann Arbor, 1943. 


p. 585, cos 29°10’55” for 87397577 read 87307577. 


W. H. RoBERTSON 
Sydney Observatory 
Sydney, Australia 


197.—U. S. Hyprocrapuic OrFice, Tables of Computed Altitude and Azi- 
muth. Publication 214. 


During the course of a recent computation project at the U. S. Naval 
Observatory it was found advisable to put on punched cards approximately 
one-third of the tabular functions contained in H.O. 214, Tables of Com- 
puted Altitude and Azimuth, volumes 1-9. The functions for latitudes 
30°-39° were taken from a manuscript, kindly furnished by the Hydro- 
graphic Office, which had been used in correcting the grosser errors in the 
first computation of volume 4. This made possible a critical examination 
of these data by differencing methods designed to detect minimum errors 
of 0/2 in altitude and 0°2 in azimuth, that is, errors of two units or more 
in the last printed place. 

Of the eighteen errata listed below only one might seem to qualify as 
a computational error, the others apparently being typographical faults. 
Even if all eighteen were considered to be computational errors, they would 
amount to only seven thousandths of one percent of the total number of 
functions examined. 

This is the most extensive study of the accuracy of H.O. 214 that has 
yet been made. Earlier conclusions [MTAC, v. 1, p. 81, v. 2, p. 182, v. 3, 
p. 139, 315] were based on relatively small samples. If the data examined 
are representative of the data not examined (a reasonable assumption) it 
may be concluded that the number of errors of two units or greater in H.O. 
214 average six per volume, which is a high standard of accuracy. 
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Errata: 
V. 1; p. 75, lat 3°, dec 13°, h.a. 82°, for alt 7°00‘2 read 7°06! 2. 
p. 131, lat 5°, dec 21°, h.a. 0°, for alt 64°30/0 read 64°00'0. 

V. 2; p. 110, lat 14°, dec 11°, h.a. 23°, for alt 67°23!1 read 67°21'3. 

. 186, lat 17°, dec 7°, h.a. 80°, for az 85°1 read 86°1. 
V. 3; p. 31, lat 21°, dec 7°, h.a. 20°, for alt 55°77!9 read 55°47!9. 

. 112, lat 24°, dec 12°, h.a. 44°, in alt 46°39/9 the 6 is poorly printed. 
V. 5; p. 55, lat 42°, dec 1°, h.a. 31°, in alt 38°41/9 the 9 is poorly printed. 


p 
p 
p 
p 
p 
p. 120, lat 44°, left hand h.a., for first 13° read 12°. 

p. 121, lat 44°, dec 30°, h.a. 31°, for az 150°0 read 153°0. 

p. 141, lat 45°, dec 17°, h.a. 14°, for alt 26°43'2 read 26°42'3. 
p 

p 

p 

p 

p 

p 

p 

p 


V. 6; p. 7, lat 50°, dec 9°, h.a. 51°, for alt 16°14!7 read 16°14'5. 

. 58, lat 52°, dec 10°, h.a. 88°, for az 95°4 read 85°4. 
V. 7; p. 41, lat 61°, left hand h.a., for second 132° read 133°. 

. 147, lat 65°, dec 30°, h.a. 146°, for alt 8°26!7 read 8°36!7. 
V. 8; p. 244, lat 79°, dec 18°, h.a. 49°, for alt 24°95/2 read 24°59'2. 
V. 9; p. 37, lat 81°, dec 18°, h.a. 160°, for alt 9°3/13 read 9°31'3. 


. 162, lat 86°, dec 11°, h.a. 62°, for alt 12°5/13 read 12°51'3. 
. 217, lat 88°, left hand h.a., for second 132° read 133°. 


Raynor L. DuNCOMBE 
U. S. Naval Observatory 


Washington 25, D. C. 


UNPUBLISHED MATHEMATICAL TABLES 


131[E, L].—R. M. CoGuian & R. C. T. Smita, Table of roots of sinz = kz. 
Typewritten Manuscript, 2 leaves, on deposit in UMT Fite and with 
Aeronautical Research Laboratories of the Department of Supply, Box 
4331 G.P.O. Melbourne, Australia. 


The table gives 6D values of the real and imaginary parts of the first 
11 zeros of sing + kz for + k = 0(.25)1. [The results for k = 1 have 
been published; see MTAC, v. 3, p. 414, RMT 611.] 


132[F].—R. A. LrENARD, Tables of the factors of 2” — n — 2 and 2" — n — 3 

for m = 1(1)1000. 

These two functions of are remarkable in that they do not seem to 
represent primes. In fact 3 is the only prime that they are known to repre- 
sent. Most values have small prime factors. In this respect these functions 
resemble CULLEN’s! function 1 + x2*. At least one prime factor of 2" — n — 2 
is given for m < 1000 except for the 17 values: m = 253, 323, 355, 455, 
493, 497, 517, 535, 559, 589, 649, 713, 749, 815, 895, 901, 979. At least one 
prime factor of 2" — nm — 3 is given for m < 1000 except for the 12 values: 
n = 162, 210, 254, 320, 330, 416, 590, 650, 738, 780, 872, 914. 


R. A. LrENaRD 
95 Rue Béchevelin 


Lyon, France 


1A. J. C. CunnincHam & H. J. Woopatt, “Factorisation of Q = (2¢+g) and 
(q-2¢ + 1),”” Messenger Math., v. 47, 1917, p. 1-38. 
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133[F].—F. L. Mixsa, Table of primitive Pythagorean triangles with their 
perimeters arranged in ascending order from 119992 to 499998. 506 type- 
written leaves on deposit in UMT Fite. 

This table is a continuation of UMT 111 [MTAC, v. 5, p. 28], a table 
by A. S. ANEMA to 120000. The introduction gives a table showing the 
number of triangles whose perimeters do not exceed P for P = 120000 
(10000) 500000. The total number of these triangles is 35114. LEHMER’s 
asymptotic formula gives 35115. Similar data are given for pairs of tri- 
angles having equal perimeters of which there are 1750. There are 65 cases 
of 3, and one case of 4 isoperimetric triangles. 


F. L. MrxKsa 
613 Spring Street 
Aurora, Illinois 


134{L].—Y. L. Luxe & D. Urrorp, Tables of fo” (st)—(2 + t —Va? + #?)e~‘#dt. 
8 mimeographed leaves on deposit in UMT FILE and also available from 
Midwest Research Institute, Kansas City 2, Missouri. 


The tables give the real and imaginary parts U + iV of the integral 
given in the title together with the function 
U + log 22 + y — 1. 
6D values of the three functions are given for z = 0(.01).1(.1)4(.2)5. 


135[L].—J. E. Wi_krns Jr. & Nina Kroporr, Table of Laguerre Functions. 
Seven mimeographed leaves on deposit in the UMT Fie. 


The table gives 4D values of 
Ly(x)/n! = M(—n, 1,2) = 5 (— 2)@)/k! 
k=0 


for m = 2(1)7 and x = 0(.1)10(.2)20. [See MTAC, v. 1, p. 361, 425, v. 2, 
p. 31, 267.] 

J. E. Wixxins, Jr. 
Nuclear Development Associates 
80 Grand Street 
White Plains, N. Y. 
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TECHNICAL DEVELOPMENTS 


Provision for Expansion in the SEAC 


In developing the SEAC, two divergent objectives had to be attained. 
The first objective was to get a modest performance high-speed computer 
into operation at the earliest possible date; the second objective was to 
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develop an effective experimental tool for evaluating advanced computer 
components and systems. In order to reconcile these conflicting aims, it 
was decided that the SEAC should be made expandable, i.e., should consist 
of a central nucleus (with a high priority on early completion) to which 
additional units could be added later, as time permitted or as initial oper- 
ating experience suggested. 

The central nucleus of the SEAC was completed and put into regular 
operation in the Spring of 1950. Since then, several additional units have 
been incorporated into the machine’s memory, input-output, and control 
equipment. 


1. Memory.—The first major change made in the SEAC was to double 
the size of the high-speed internal memory. Initially the SEAC was oper- 
ated with an acoustic delay-line memory holding 512 words stored serially 
in 64 mercury tanks (with average access time: 168 microseconds). Subse- 
quently an additional high-speed memory was incorporated which holds 
512 words stored in 45 Williams’ tubes. All 45 tubes simultaneously supply 
the computer with 45 binary digits in parallel (with average access time: 
12 microseconds). Although the SEAC is in other respects a serial machine, 
provision was made in the original design for the use of both serial and 
parallel types of memory. Conversion of the word-form back and forth 
from serial to parallel is accomplished by means of a shift register with a 
special flexible control. 

2. Input-Output.—The same shift register equipment is sufficiently ver- 
satile to allow the SEAC to operate with completely asynchronous, serial 
or parallel, input-output equipment. For example, the SEAC has been 
operated in conjunction with Teletype punched-paper tape, with Teletype 
keyboard and printer, with a visual (cathode-ray tube) high-speed graph 
plotter, with single-channel magnetic wire, and with single-channel mag- 
netic tapes. Additional units which are planned to be incorporated in the 
near future include multi-channel magnetic tapes, Flexowriter punched- 
paper tape and printing equipment, punched-card-to-magnetic tape con- 
version equipment, and a high-speed photographic printer. Also, an external 
selector is being added to the SEAC which will be capable of selecting (under 
the computer’s control) from among any of ten different input-output units. 

3. Control—Another major change to be made in the SEAC will be 
the addition of equipment to the control. Provision was made for this 
expansion in the original design of the machine, and, as a result, the operator 
of the computer will be able to choose (before entering his problem) between 
two distinct modes of operation (i.e., between the 4-address or the 3-address 
mode) merely by throwing a few toggle switches. 


The main distinction between the 4-address and the 3-address mode is 
indicated in Table I, which lists the information contained in an instruction- 
word in each system. Note that the first three items of information are 
the same in both systems. These items indicate where in the memory the 
operands are located (a and 8), and also where in the memory the result 
of a given operation on them is to be sent (7). In the 4-address system, the 
instruction also contains an explicit statement on where the next instruction 
is located in the memory (4). In the 3-address system, however, there is no 
such explicit statement. Instead, instructions in this system are automati- 
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cally sequenced, i.e., successive instructions being located, ordinarily, in 
successively-numbered memory positions. 

The types of operations which can be performed in both the 3-address 
and 4-address systems, are similar. They include (besides input-output 
cperations) addition, subtraction, multiplication, division, and a type of 
extract operation called “logical transfer.’’ Also, three types of comparison 
operations (sometimes called discrimination or branch operations) can be 
performed. In a typical comparison operation, when the result of the com- 
parison between two operands indicates that the first operand is greater 
than or equal to the second operand, then the next instruction is chosen 
normally, in the manner just described. But, if the second operand is the 
greater, then the next instruction in both systems is chosen from address y+. 

4. Floating Address Feature.—In the 3-address system, a special auto- 
matic feature is included which is believed to be novel and which is referred 
to as the Floating Address feature. This was incorporated into SEAC in 
order to shorten and simplify the work of the programmer. It does not 
accomplish any effects which the human programmer could not also accom- 
plish at the cost of additional labor. The scheme does, however, automatize 
certain time-consuming programming operations, thereby completely relieving 
the programmer of any need for concerning himself with them. This is ac- 
complished at the cost of a relatively small amount of additional equipment. 

In brief, this feature makes it easy for the programmer to reutilize 
standard routines or lists of instructions which have already been compiled 
for carrying out certain frequently-used operations, such as taking a square 
root, making a binary-to-decimal conversion, and so on. Obviously he 
would like to make full use of any such material originally prepared for 
one problem in the solution of other problems. But there is a difficulty in 
the way of such straightforward re-use of standard routines. This arises 
from the fact that the memory locations into which he must insert the stand- 
ard list of instructions will be different in general for each new problem. 
This difference forces him to change the memory address by which he 
identifies each individual instruction in the routine each time he uses the 
routine in a different problem. The floating address feature performs this 
change of address automatically each time the routine is used, and in a 
manner which requires no attention on the part of the programmer. 

Table II illustrates the characteristics of the floating address feature. 
In preparing a routine originally, the programmer can designate any address 
as either absolute or relative. An absolute address is interpreted in the usual 
sense merely as a number identifying a specific memory location. A relative 
address, however, identifies a position in the memory by specifying its 
displacement relative to the position of the instruction-word itself. To 
illustrate how this works: if an instruction located in memory position m 
specifies that the address of the first operand is number 17 relative, then the 
control will take the word in memory location m + 17 and use it as the first 
operand. To refer to location m — 17, the programmer uses as address the 
complement of 17 (2 — 17, relative). It is possible to use any combination 
of absolute or relative addresses in any instruction word. 

By utilizing this feature, therefore (by using relative-address designa- 
tions for all references to instructions inside the routine itself), routines 
may be prepared in advance once and for all, transcribed on tape, and 
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stored away in the library files. For any subsequent problem, then, each 
routine can be inserted in the memory in any arbitrary location whatever, 
and executed immediately without need for any preparatory modification to 
adapt it to the particular memory location in which it happens to fall. 

The second major characteristic of the floating address feature is intended 
to relieve the programmer of another related type of burden. To illustrate 
this, let us represent the various storage locations of the internal memory 
by the boxes shown on the flow chart (Table III). Each box contains an 
instruction. Each group of boxes represents a typical routine, a list of 
instructions stored in consecutively numbered address locations. The right- 
hand list represents a sub-routine describing some standard process which 
has been entered into the memory verbatim from a strip of library tape. 
The left-hand list represents the main routine prepared specially for the 
problem at hand. During the course of any problem, the control usually 
follows the main routine for a while, then jumps to some sub-routine, then 
jumps back to the main routine, and so forth. Now the same sub-routine 
generally will be referred to repeatedly during the course of a problem. 
(E.g., square roots may be sprinkled in many places throughout the com- 
putation.) The point of departure from the main routine will in general be 
different in each case. Under the 4-address SEAC system, the programmer 
has to arrange, before embarking on the sub-routine, to plant some sort of 
clue in some pre-arranged location (or at the end of the sub-routine itself) 
which will lead the control back to the main routine at the proper point of 
return each time. This means programming extra operations. 

With the 3-address floating address system, the return to the proper 
point in the main routine always takes place automatically and requires no 
attention at all from the programmer. This is accomplished in the following 
manner: The system provides two small storage registers, each of which is 
capable of holding an address number. Normally, as soon as any instruction 
is executed, the contents of one of these registers are advanced one unit, and 
the number it contains is taken as the address of the next instruction. 
The contents of the other register are held unchanged. A single binary digit 
(called “‘d” on Table II) is reserved in each instruction word for designating 
which one of the two registers is to be referred to. 

To be more specific, on main routines the instructions are normally 
written with the register designation digit d = 0. These main routines are 
sequenced by the counter register called #0. To embark on a sub-routine, 
the programmer writes a special “‘jump’’ instruction (with d = 1) which 
causes the other counter register (called #1) to take over. The sub-routine, 
always written with the digit d = 1, continues using counter #1. The final 
instruction of the sub-routine, however, is written with digit d = 0. 
Immediately after executing the final instruction, therefore, counter #0 is 
referred to for the address of the next instruction, and the control is back 
on the main routine at the point immediately following where it left off. 
All sub-routines are written in an invariant form in this respect, and no 
other attention on the part of the programmer is required. 

Some of the alternative courses followed by the control are indicated in 
Table III: normal operation, jump operation, and the path of automatic 
return. By combining this dual-counter feature with the relative address 
feature, a number of other useful operations may be performed. As an 
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example of one: by a jump operation with a relative resetting for counter #0 
at the end of the sub-routine, automatic return may be made to take place 
to any part of the main routine. 

Now the question naturally arises: How much does this added control 
equipment increase the size of the computer? The answer in terms of the 
total tube count of the complete computer is: about 3% (30 tubes out of a 
total of about 1050). Since the SEAC was originally designed with the 
possibility of these additions in mind, the actual process of attaching new 
units is quite simple and straightforward. 

Is annexation of such additional equipment justified? We believe it is. 
The value of the computer ought to be measured by its total productivity 
in terms of complete problems solved. Experience to date with the SEAC 
suggests that in the long run the main bottleneck on problem output may 
be the ability of the human operator to prepare and check adequately the 
programs that he enters into the machine. Preparing a program for a large 
scale problem is a highly intricate process and demands complete accuracy 
of detail. Furthermore, checking a program which appears to be going 
astray on the machine often turns out to be very time-consuming. It is 
found, for example, that the time necessary for locating and correcting 
program errors is often comparable with the time spent in locating and 
correcting machine errors. With the advent of (first) more complicated 
problems and (second) more reliable machines, the balance will become 
even more heavily weighted against the programmer. By simplifying his 
job and reducing the number of details with which he has to concern himself, 
we think we can reduce the incidence of human error and thereby gain 
markedly in efficiency. For this reason, we feel that the supplementary 
equipment just described will contribute substantially towards increasing 
the SEAC’s problem-solving productivity. 


NBS Electronic Computers Lab. 
Washington, D. C. 


ALAN L. LEINER 


TABLE I. THE SEAC’S CONTROL 


NOW IN OPERATION BEING INCORPORATED 
4-address system 3-address system 
(with Floating Address feature) 
Each instruction word contains: Each instruction word contains: 
(1) Address of first operand (1) Address of first operand 
(a) ... 10 binary digits ; (a) ... 12 binary digits 
(2) Address of second operand (2) Address of second operand 
(8) ... 10 binary digits (8) ... 12 binary digits 
(3) Address of result (3) Address of result 
(y) ... 10 binary digits (y) ... 12 binary digits 
(4) Address of next instruction (4) Floating address information 
(6) ... 10 binary digits @,b,¢.d) ..... 4 binary digits 
(5) Symbol specifying type of operation (5) Symbol specifying type of operation 
to be performed to be performed 
.... 4 binary digits .... 4 binary digits 
(6) Halt signal (6) Halt signal 
.... 1 binary digit .... 1 binary digit 
Total: 45 binary digits Total: 45 binary digits 
in each word in each word 








DU 
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TABLE II. FLOATING ADDRESS INFORMATION 


RELATIVE ADDRESS INFORMATION 
(a) indicates whether a is an absolute or relative address .... 


1 binary digit 

(b) indicates whether 6 is an absolute or relative address .... 
1 binary digit 

(c) indicates whether 7 is an absolute or relative address .... 
1 binary digit 

DUAL REGISTER INFORMATION 
(d) modifies choice of next instruction ......... 1 binary digit 
THR, vkcas 4 binary digits 


TABLE III. CHOICE OF NEXT INSTRUCTION IN THE 
THREE-ADDRESS (FLOATING) SYSTEM 


POSSIBILITIES JUMP operation, indicated by: — — — > 
NORMAL operation, indicated by: ——— > 
PATH OF AUTOMATIC RETURN, indicated by: ==> 
















































































MAIN ROUTINE SUB-ROUTINE 
Address No. Address No. 
(er 105 Branch Instruction 7 7 368 Initial Instruction 
C / 
106 / 369 
/ 
107 / 370 
/ 
108 4 371 Branch Instruction 
4 
109 Jump Instruction F~ 372 ) 
110 373 
111 374 
YG 112 ) 375 Jump Instruction 
ar 
113 et 
rg 
” 
-~ 
114 ” 
115 
Etc. 











Notes on Numerical Analysis—6 
The Convergence of Seidel Iterants of Nearly Symmetric Matrices 


If A is a matrix, we may write A = D + L + R. Here D is the matrix 
of its leading diagonal elements, L the triangular matrix of the elements 
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below this diagonal, and R the triangular matrix of the elements above. 
If B is the Seidel! iterant of A, then B = — (D+ L)“R. 
It is known that if A is symmetric, lim B* = 0 if A is positive definite.” 


n> 0 

By continuity one may infer that if A is ‘nearly equal’ to some symmetric 
matrix C, the Seidel iterants of A and C would behave in a similar way 
with regard to convergence. For, although there is a limiting process in- 
volved, the convergence of B* to 0 depends on the maximum modulus of 
the characteristic roots of B, and this is a continuous function of the ele- 
ments of B. The problem attempted in this paper is to determine a precise 
meaning for the notion ‘nearly equal.’ : 

Any matrix A can be written in the form A = C + S, where C is sym- 
metric and S = \ — 0’ is skew-symmetric (d is defined in the Theorem). 
C may be regarded as the symmetric matrix nearest to A, and \ is then a 
measure of its asymmetry. Our solution of what is to be understood by 
‘nearly equal’ is given in terms of C and X. The result is simplest when 
D =I, the unit matrix. We state the result in that case as a theorem. 
The minor adjustments required when D = J, all a;; > 0, are explained at 
the end of the paper. 

Theorem. If A=I+L+R, C=1I+}(L+4+ R’) + 34(L' + R), 
B= —-—-([+L)7°R, }\=}(L—R’), 


M=I+dA4+NV—-F{A-M)T+L4+R)77+L—-R) 
+(U4+L'’-—-R)I4+ L' + RF’) - \)}, 
then, provided M is positive definite, lim B* = 0, or # 0 according as C is 


or is not positive definite. 
For any matrix A, let [N(A) }? = 2;,0%;. If we write M =I + M*, 
then it may be proved that M is certainly positive definite if N(M*) < 1. 
From known inequalities for the norm of a matrix,’ 
N(M*) < 2N(A) {1+ NCI + 24+ R)“)INT + 2 + R)}. 


Hence if an estimate of N{(J + L + R)} may readily be obtained, the 
theorem gives a practical test for the convergence to 0 of the Seidel iterants 
of a nearly symmetric matrix. 
We proceed with the proof of the theorem. Let x be an arbitrary 
vector ; let 
x= Bx, y= Br, s=x-y. 


We suppose J — B to be non-singular, and so obtain 
x= (J — B) "sz, «+ y= (7+ B)U — B)z 


= (I — B)"(I + B)z. 
Since C is symmetric, ' 


(1) x'Cx — y'Cy = (x’ — y')C(x + y) = oC — B)“( + Biz. 


Since B = — (I+ L)“R, and C=A+2X —), substituting in (1) 
we get 


(2) x/Cx-—yCy=27(7+L+R-A4t+N)I4+L+R704+L — Riz. 


We interrupt the proof to state a lemma: 








lin 


n—- 


Sin 
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Lemma. If C is a symmetric matrix, y = Bx, x'Cx — y'Cy > 0, then 
lim B* = 0 or ¥ 0, according as C is positive definite or not. 


A formal proof of the lemma is given elsewhere.* 
If now A = 0 then (2) gives 


x’Cx — y'Cy = 2s>0 if e #0. 


Since z = (I — B)x and I — B is not singular, z + 0 if x ~ 0. By the 
lemma the known result for the convergence of the Seidel iterants then 
follows immediately. 

If A is not symmetric, we may continue proving the theorem as follows. 
For any matrix P, 2’Pz = $2/(P + P’)z, and 3(P + P’) is a symmetric 
matrix. Applying this result to the matrix on the right of (2) we get 
x’Cx — y’Cy = 2’ Mz, where M is the matrix defined in the enunciation of 
the theorem. The theorem now follows from the lemma. 

IfA = D+L+R, and D # I, all a, > 0, the theorem is still true if 
we generalize somewhat the definitions of C and M. Let 6 be the diagonal 
matrix whose elements are the positive square roots of the elements of D. 
Then & = D. The Seidel matrix of A is then 


B= — (D+ LOR = — (84+ DLR = —- H+ HLH) FR. 


If we write B* = 6Bé5-, the characteristic roots of B* are those of B. 
Hence the B iterants will converge to zero or not according as the B* 
iterants converge to zero or not. Again B* = — (J + 6“*Lé")“"d7 RO. 
Hence B* is the Seidel iteration matrix for J + 6“ Ld“ + 6° RO. 

Now, if 7+ 2+R is symmetric, so is J + 6"Lé" + 5“R5-.. The 
theorem may thus be extended to the more general case provided that in 
the definitions of C, B, \ and M, we replace L and R by 6“ Lé— and 6" Rd“ 
respectively. 

If A=D+L+R is a symmetric matrix, the condition that A is 
positive definite in order that the Seidel iterants of its associated matrix 
converge to zero is not only sufficient but, in the case when the elements of 
D are all positive, is also necessary.® For we observe that in the case D = J 
this result is a special case of the theorem. If D ¥ I, we form the matrix 
A* = 5"As— where 6 is defined as above. If now M is any leading minor 
of A and M* the corresponding minor of A*, and 4,, the elements in 6 corre- 
sponding to M, then M* = 6,,-'M6,,— and so det (M*) = (det 4,,~")* det (M). 

If the elements in D are all positive, (det 4,,~')? is positive and det (M*) 
and det (M) are both of the same sign. Hence A and A* are both positive 
definite, or neither is positive definite. Since the Seidel iterants of the 
matrices corresponding to A and A* either both converge to zero or neither 
does, and the leading diagonal in A* is J, the result for A follows from what 
has been proved for A*. 

I wish to thank OLGa Taussky-Tonp for the encouragement I received 
from her in the course of the preparation of this paper. I also wish to thank 
my colleague, S. E. Cruise, for reading and checking the manuscript. 


P. STEIN 
University of Natal 


Durban, Union of South Africa 
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1 See, for example, HiLpa GErRINGER, “On the solution of systems of linear equations 
by certain iteration methods,” Reissner Anniversary Volume, Contributions to Applied 
Mechanics, Ann Arbor, 1948, p. 365-393. 

2 This is Seidel’s result. For a recent proof see R. von Mises & HILDA GEIRINGER, 
“Praktische Verfahren der Gleichungsauflésung,”’ Zeitschr. f. angew. Math. u. Mech., v. 9, 
1929, p. 58-77, 152-164. 

3 See, for example, p 1044 in Joun von NEuMANN & H. H. GotpstineE, “Numerical 
a of matrices of ‘high order,”” Amer. Math. Soc., Buli., v. 53, 1947, p. 1021-1099, 
TEIN, ‘“‘Some general theorems on iterates,’’ to appear in NBS Jn. “of . Res. 

bd A different proof of this result is given by E. REICH, * ‘On the convergence of the classi- 
cal —— of solving linear simultaneous equations,’’ Ann. Math. Stat., v. 20, 1949, p. 
448-451 


BIBLIOGRAPHY Z-XVII 


1. Anon., ‘“‘Conversion between analogue and digital data,’””» ONR London, 

Technical Report ONRL-73-50, July 1950, 4 pages. 

This brief report from the Scientific Liaison Officer of the American 
Embassy in London, England, summarizes a symposium held at the Mili- 
tary College of Science, Shrivenham, 5-6 July, 1950. The following pre- 
pared papers were presented and discussed: 


(1) “A method of telerotation by pulse code modulation,” by R. H. 
BARKER of the Signals Research and Development Establishment, 
Christchurch, Hants., 

(2) “The production and use of binary coded discs for data trans- 
mission,”’ by W. S. Etttott, R. C. Rossins, and D. S. Evans of 
the Research Laboratories of Elliott Bros., Borehamwood, 

(3) ‘‘Telerotation by parallel binary codes,’’ by E. A. JoHNsoN of the 
Radar Research and Development Establishment, Great Malvern, 
Worcs., 

(4) “A servo system operating on discontinuous information,” by C. 
Ho.t-SmitH and D. LAWDEN of the Faculty of Instrument Tech- 
nology, Military College of Science, Shrivenham, 

(5) ‘‘Note on digital computing for the extrapolation of discontinuous 
information,” by G. C. Toorit of the Faculty of Instrument Tech- 
nology, Military College of Science, Shrivenham, 

(6) “Some methods of conversion of data from analogue (shaft rotation) 
to digital form,’’ by K. V. Drprose of the Royal Aircraft Establish- 
ment, Farnborough, Hants., 

(7) ‘The electromagnetic cam: a technique for wave manipulation using 
the nonlinear properties of magnetic materials,’’ by W. LAWRENCE 
of the Signals Research and Development Establishment, Christ- 
church, Hants. 


Copies of these papers have been forwarded to Code 250, Office of Naval 
Research, Washington 25, D. C., and the proceedings of the Symposium 
(except for these papers) will be prepared and distributed by the Ministry 
of Supply. 

The different systems described for conversion between shaft rotation 
and digital codes had different chief features such as compactness, accuracy 
of one part in 20,000, or a three microsecond digit time. The prediction 
problem was treated both by direct use of the discrete data and by poly- 
nomial approximation to the prediction function. 


R. D. ELBourN 
NBSCML 
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2. ANON., ‘Desk-sized electronic computer developed by Northrop,” Ma- 
chinery, v. 57, Mar. 1951, p. 168. 


Short expository article on the MADDIDA. 


3. ANON., Digital Computer Newsletter, ONR Mathematical Sciences Divi- 
sion, v. 3, April 1951, 6 pages. 
The contents are as follows: 


1. The ERA 1101 Computer 

2. Moore School Automatic Computer 

3. The Institute for Advanced Study Computer 

4. Aberdeen Proving Ground Computers: The ENIAC, The 
EDVAC, The ORDVAC 

. National Bureau of Standards Western Automatic Computer 
(SWAC) 

. National Bureau of Standards Eastern Automatic Computer 
(SEAC) 

. Project Whirlwind 

. Relay Digital Computer, Imperial College, Univ. of London 

. Zuse Computer Model IV, at Zurich, Switzerland 

. Data Handling and Conversion Equipment: Stavid Engineering 
Data Conversion Equipment, Signal Corps Angular Position 
Encoders, and Zatacoding 
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4. ANON., Digital Computer Newsletter, ONR Mathematical Sciences Divi- 
sion, v. 3, July 1951, 6 pages. 
The contents are as follows: 


General Purpose Computers 


1. Burroughs Laboratory Computer 

2. UNIVAC Dedication 

3. The Raytheon Digital Computer 

4. The National Bureau of Standards Eastern Automatic Computer 

5. Whirlwind I (M.I.T.) 

6. The ORDVAC (Univ. of III.) 

7. Moore School Automatic Computer (MSAC) 

8. The Institute for Advanced Study Computer 
Special Purpose Computers 

1. The CRC 101 Digital Differential Analyzer 

2. The MADDIDA 


Auxiliary and Conversion Equipment 


1. Flying Typewriter 
2. Conversion and Display Equipment 


5. Anon., ‘““Guided-missile computer,” Electronics, v. 24, Apr. 1951, p. 138. 
This is an expository paper which describes the hybrid analog-digital 


computer being developed at Project Typhoon, Special Devices Center of 
ONR, Port Washington, N. Y. 
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6. T. J. Conno.zy, S. P. FRANKEL, & B. H. Saaz, “Application of auto- 
matic digital computing methods to the production of phase behavior,” 
Electrical Engineering, v. 70, Jan. 1951, p. 47. 

Punched card methods using the type 604 calculating punch were devel- 
oped to predict the thermodynamic properties of homogeneous and hetero- 
geneous multicomponent hydrocarbon systems by solving the Benedict 
equation of state. Methods described permit evaluating the composition 
and the molal values of volume, entropy, and enthalpy of the coexisting 
phases at one heterogeneous state in about one hour of computing, or for 
a single homogeneous phase in about ten minutes. Thus computing par- 
ticular values as required is feasible whereas tabulating completely the 
thermodynamic properties of multicomponent heterogeneous systems is not. 
A bibliography of 36 titles is included. 


R. D. ELBourn 
NBSCML 


7. D. H. GripLey & B. L. SarRAHAN, “Design of the Naval Research Lab- 

oratory Computer,” Electrical Engineering, v. 70, Feb. 1951, p. 111. 

It is expected that the Naval Research Laboratory Electronic Digital 
Computer (NAREC) will be working by January 1952 and will be in pro- 
ductive operation within six months following. It will have 1024 words of 
electrostatic storage and 2048 of magnetic drum storage. Words will con- 
sist of 44 binary digits and sign; instruction words will be interpreted as 
pairs of one-address orders. Input and output will be on magnetic tape. 
The machine will have a multiplication time of 300 microseconds. 


MICHAEL MONTALBANO 
NBSCL 


8. LAWRENCE P. LEssING, ‘‘The electronics era,” Fortune, v. 44, July 1951, 
p. 79-83, 132-138. 


Briefly mentioned in this article are the important recent developments 
' in the computer field. To quote from the paper: “It may be that the most 
significant accomplishment for electronics for the future is its ability to 
count.” 


9. D. M. McCatium & J. B. Situ, ‘Mechanized reasoning, logical com- 
puters and their design,’’ Electronic Engineering, v. 23, Apr. 1951, p. 
126-133. 

A small electrical computer, built at Ferranti, Ltd. in Edinburgh, to 
evaluate Boolean propositional functions of up to seven binary-valued 
variables is described. A function is set up by plugging together logical 
connective boxes which contain relays, selenium diodes, and indicator lamps. 
The computer then tries automatically in sequence the 128 combinations 
of values of the variables and stops on any combination for which the 
function is ‘‘true.”’ 

One possible extension described is particularly interesting because a 
computer using it might find reasonably quickly a solution for a function 
of so many variables that systematic scanning would take a prohibitively 
long time; moreover, the process appears to be more like that used by the 
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human brain. The method is to assume some set of values of the variables, 
which usually will not be a solution, then to identify those variables whose 
values appear to be responsible for the false result and, by feedback, to 
cause them to change more or less at random until a true result is achieved. 
A scheme for avoiding circular sequences of changes is suggested. 


NBSCML R. D. ELBourn 


10. AmBros P. SpPEIsER, Entwurf eines Elektronischen Rechengerates, Inst. f. 
angew. Math., Mitt., Ziirich, No. 1, 1950. 


This is a description of the logical design of a digital computer intended 
to be built at low cost and to operate at relatively low speed, but otherwise 
to equal the performance of most of the modern automatic machines. No 
new engineering development is contemplated; use is made only of existing 
and well-tested components, especially of those developed by H. H. Aiken. 
The memory organ is a magnetic drum holding 1200 words, each consisting 
of sign and 12 decimal digits. The average access time is 16 milliseconds. 
The arithmetic organ operates serially, with a pulse repetition rate of 40 k.c. 
The addition time is .8 ms; the average multiplication time 30 ms. Multi- 
plication is accomplished by repeated addition, division by an original 
modification of the repeated subtraction process. Instructions are written 
in a single-address system, pairs of addresses being stored in the memory 
as single words. The arithmetic unit and part of the control are electronic, 
but address selection is performed by relays. The design calls for 1000 
vacuum tubes and 500 relays in the entire machine. Teletype equipment is 
used for input and output, without provision for external storage. There 
is no automatic checking, except that the machine stops when certain 
“forbidden” pulse combinations occur. 

NBSCL FRANZ L. ALT 


11. H. W. Spence, “Systematization of tube surveillance in large scale 
computers,” Electrical Engineering, v. 70, July 1951, p. 605-608. 
Ninety per cent of the service interruptions of ENIAC at Aberdeen 

Proving Ground are caused by vacuum-tube failures. This paper describes 

the first five months’ operation of a new tube-surveillance program. The 

average life of some 5000 tubes discarded varied among ten types from 

5,500 to 12,500 hours. Much of the difference between types is attributed 

to different operating and testing conditions. For most types low emission 

is the predominant cause of failure, especially so in normally OFF tubes. 

If only low-emission failures are considered, the life of tubes used in the 

ON condition is better than that of the same type used in the OFF condi- 

tion by 50 per cent or more. Poor cutoff, short circuited elements, and 

burned out filaments may be considered accidental failures because life ter- 
minated by these averaged less than half as long as that for low-emission 
failure. 

Because it contains the first extensive tube life data appropriate to 
computer design, this report is especially welcome and further results are 
eagerly awaited. 

NBSCML R. D. ELBourN 
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NEws 


Department of the Air Force and National Bureau of Standards.—A symposium on 
Linear Inequalities and Programming was held in Washington, D. C., on June 14, 15, 
and 16, 1951. This meeting was sponsored by the Department of the Air Force and the 
National Bureau of Standards. The program was as follows: 


Thursday, June 14, Morning 

Linear inequalities 

Welcoming talk 

Theorems of alternatives for pairs of mat- 
rices 

Remarks on the history of work on linear 
inequalities 

Optimization in unitary spaces 


Polyhedral cones 
Geometrical significance of the simplex 
method 
Discussion 
Thursday, June 14, Afternoon 
Applications 


Choice models and non-choice models 

Research program of Project SCOOP 

Discussion 

Friday, June 15, Morning 

Computational theory and techniques 

Welcoming talk 

New techniques for linear inequalities and 
optimization 

Methods of solving linear equations 

Gradient methods in Lagrangian prob- 
lems and their game theoretical inter- 
pretation 

Properties of Leontief matrices 

Problems of formulation 


Efficiency prices for decentralized deci- 
sions 

Suitability scales for allocation problems 

Personnel classification 


Design of an optimal battery of tests 
Saturday, June 16, Morning 

Nonlinear and linear problems 

Theory of the transportation problem 

Convex programming 

Least ballast shipping required to meet a 
specified shipping program 

Problem of contract awards 

A gasoline blending problem 

Discussion 


J. Topp, Chairman, USAF 
Gen. F. J. Dau, USAF 
A. W. Tucker, Princeton Univ. 


T. S. Morzxin, NBSINA 


E. W. BaRANKIN, Univ. of Calif. and RAND 
Corp. 

A. Cuarnes, Carnegie Inst. of Tech. 

J. P. MayBerry, Princeton Univ. 


T. C. Koopmans, Chairman, Cowles Com- 
mission, Univ. of Chicago 

W. W. Leontier, Harvard Univ. 

M. K. Woop, USAF 


G. B. DantziGc, Chairman, NBSCL 
Gen. F. J. Dau, USAF 
T. S. Morzxin, NBSINA 


G. ForsytHE, NBSINA 
L. Hurwicz, Univ. of Illinois, and Cowles 
Commission, Univ. of Chicago 


Y. K. Wons, Princeton Univ. 

P. A. SAMUELSON, Chairman, MIT and 
Ranp Corp. 

T. C. Koopmans and G. DEsBrREv, Cowles 
Commission, Univ. of Chicago 

J. L. Hottey, USAF 

D. F. Voraw, Jr., Yale Univ., and A. 
OrpEN, USAF 

M. A. Woopsury, Princeton Univ. 


E. D. SCHELL, Chairman, USAF 

M. M. FLoop, Ranp Corp. 

D. W. BLAcKETT, Princeton Univ. 

I. HELLER, Logistics Research, George 
Washington Univ. 

L. Go_pstern, USAF 

W. Cooper, Carnegie Inst. of Tech. 
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Saturday, June 16, Afternoon 


Computational theory and techniques F. L. Att, Chairman, NBSCL 
Matrix inversion H. H. Gotpst1n_, Inst. for Advanced Study, 
Princeton Univ. 
A unified technique for matrix inversion, A. OrpEN, USAF 
linear inequalities, games, optimization 


A short proof of the dual theorem G. B. Dantzic, USAF 

The projective method and computational C.B. Tompkins, Logistics Research, George 
experience Washington Univ. 

Discussion 


Eckert-Mauchly Division, Remington Rand Inc.—The first UNIVAC, which was 
constructed for the Bureau of the Census, under NBS contract with Remington Rand, 
was delivered by Epwarp U. Connon, Director of the NBS, to Roy V. Peet, Director 
of the Bureau of the Census, on June 14th. 

During the elaborate ceremonies which took place in Philadelphia on that occasion, 
the UNIVAC was busy on a typical task of the Bureau of the Census—the tabulation of 
information about JoHN Q. CiT1zEN. In one-sixth of a second, the machine classified an 
individual according to place of residence, sex, color, age group, nationality, citizenship, 
extent of education, educational activity, farming activity, government activity, marital 
status, domestic status, veteran status, employment status, occupational group, industrial 
group, earning capacity, etc. For each county, the resulting tabulations were printed out 
by the UNITYPER under the appropriate headings. Any interested guest at the dedication 
ceremonies had the opportunity to check the crossfootings testing the accuracy of the 
machine’s operations, 

Even at the amazingly rapid rate of six persons per second, the classification of 150,000,000 
individuals and the preparation of the various and sundry tabulations is too formidable a 
task for a single machine. It is planned therefore to continue the use of punched card 
equipment as heretofore, in conjunction with the tabulation by the UNIVAC of the resi- 
dents of a few selected states. The Bureau of the Census estimates that about ten per cent 
of the total task will be accomplished by electronic equipment by the time the 1950 
tabulations will have been completed. 

During the morning of June 14th, a press seminar was held at which the following 
program was presented: 


Welcome ALBERT GREENFIELD, President, Chamber 
of Commerce, Philadelphia 
Description of the UNIVAC J. Presper Eckert, Eckert-Mauchly Divi- 


sion, Remington Rand, Inc. 
Description of Demonstration Problem James L. McPuHeErson, Bureau of the Census 
Demonstration of the UNIVAC S. N. ALEXANDER, NBS 
Press Forum 


During the afternoon the dedication ceremonies were held. The program was as follows: 


Chairman ALBERT GREENFIELD, President, Chamber 
of Commerce, Philadelphia 

Economic Implications James H. Ranp, President, Remington 
Rand, Inc. 

Official Delivery of the UNIVAC Epwarp U. Connon, Director, NBS, and 
Roy V. Pret, Director, Bureau of the 
Census 


Dedication CHARLES SAWYER, Secretary of Commerce 
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OTHER AIDS TO COMPUTATION 


A Graphic Solution of Simultaneous Linear 
Equations 


This note presents a graphic method of solving sets of simultaneous 
linear equations which can easily be extended to the inversion of matrices 
and the calculation of determinants. Its main advantages are speed and 
simplicity, its main disadvantage is the approximate nature of graphic 
solutions. 

The procedure to be followed involves the repeated use of the following 
diagram: 
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Figure 1 


The dotted lines are drawn parallel, the vertical distance between the 
top two being h;, between the bottom two being h,. Since c is parallel to 
a and since ¢ is parallel to b, the ratios c/a and e/b are each proportional to 


he/ (hy + he). That is, 
c/a = In/(y + In) = e/b. 
Therefore, e = cb/a, and 
(1) f=d-—e=d-— c/a. 
Now, if the given equations are 

OX + AyxX2 + +++ + Ainkn = Ay 

AayX1 + GaeX2 + +++ + AonXn = An 
(2) F eda geapanngy edtenlad ad 
OniX1 + AnoX2 + +++ + AnnXn = Ano 


and if Figure 1 is so drawn that the lengths of lines a and b measure ay 
and a,;, and the lengths of lines c and d measure a, and a,;, we have 


(3) f = diya = ay — A1j0i3/ aun 


where 4;;,1) is the coefficient of x; in the i-th reduced equation obtained by 
the elimination of x:. 

The graphic procedure followed here is exactly parallel to the arithmetic 
of the Gauss method and its numerous variations such as the Doolittle 
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method. Equations (2) are reduced to 


Qe2(1)X2 + +++ + denayXn = Aq 
(4) GsoyX2 + +** + AanayXn = Ana) 


OnoyX2 + +++ + AnngyXn = Snoq)- 
The process, after  — 1 repetitions, gives 
(5) Aan(n—1) Xn = Ond(n—-1)s 


which is easily solved for x,. 

Finally a “‘back solution’’ gives in turn the values x,_1, Xn-2, **-, X1- 

This entire process can be carried out graphically by repeated use of 
diagrams like that in Figure 1. We shall illustrate the process by the follow- 
ing example. 

Consider the following set of two linear equations—typical of a problem 
of regression in three variables; 


(6) Biz.s + 0.55 Bis.2 = 0.82 
0.55 Biz.s + Bis.2 = 0.68 


The solution of this problem is given in Figure 2. 
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Figure 2 


The coefficients in the equations are plotted as heavy horizontal lines 
to some convenient scale. That is, the lengths of the lines in the first row 
are ab = ay = 1; cd = dy = 0.55; and ef = ay = 0.82. The lines in the 
second row represent gh = dy = 0.55; ij = dm = 1; and kv = do = 0.68. 
The dotted lines are used to evaluate the reduced coefficients. First, we 
draw line bh and extend it to meet the O-line at m. Second, we draw mo 
perpendicular to the 0-line. Third, we draw dj and extend it to meet mo 
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at n. Fourth, we draw nc, intersecting 7j at p. Fifth, we draw fv and extend 
it to intersect mo at o. Sixth, we draw eo, intersecting kv at g. Now, 
G22) = tp, dasa) = kg, and Bis.2 = kq/ip. 

We can evaluate §13.2 graphically by drawing the heavy dotted lines 
below the segments representing the second equation. These lines are 
located on the 0-line by the intersection of line cu and a line parallel to the 
0-line at a distance 1. The value of 13.2 is measured by line st. 

A “back solution” is shown also for By.3. First, draw line rd and extend 
it to meet the 0-line at u. Second, draw uu’ horizontal. Third, draw line 
tf and extend it to intersect uu’ at u’. Fourth, draw su’, intersecting ef at w. 
The value of Bis.3 is ew. 

Whether the computations are carried out arithmetically or graphically, 
the number of steps in an n-variable problem increases rapidly as m increases 
about in proportion to m*. A diagram for a problem of four or five variables 
is much larger than Figure 2 and many more points must be located. But, 
as in the Doolittle method, the entire process is a repetition of the steps 
taken in Figure 1. And after a little practice it is unnecessary to draw all 
the lines indicated in Figures 1 and 2. All we need is to locate and mark 
points of intersection, using a transparent straightedge. 

We have found that the reduction of each coefficient takes about five 
seconds if there are no complications. The complications that arise are 
caused by a lack of space if two of the initial coefficients are almost equal, 
causing near parallelism of two lines which must intersect. This can easily 
be rectified by first going through one of the equations involved and halving 
each of its coefficients and the constant term. 

The scale should be chosen so that there is no length of as long as 8 
inches. Generally the result will have accuracy of three significant figures, 
occasionally less. It should be noticed in the computation that negative 
coefficients are measured to the left of the 0-line, positive to the right. 
If more accuracy is needed this process gives an excellent first approxima- 
tion to use in one of the iterative methods discussed by HOTELLING.! 

The equipment necessary includes a good drawing board, a T-square 
and a triangle, together with a rule, preferably transparent with several 
scales. If one uses colored pencils for consecutive reductions of the set of 
equations one can do at least a5 X 5 set without tracing the partial solution 
to another sheet. 

Joun E. MAXFIELD 
University of Oregon FREDERICK V. WAUGH 
Eugene, Oregon 

1 HaROLpD HOoTELLING, “Some new methods in matrix calculation,” Annals Math. Stat., 

v. 14, 1943, p. 1-34. 


BIBLIOGRAPHY Z-XVII 


12. H. T. Avery, “The design of office machines—such as calculators,” 
Mech. Engineering, v. 71, 1949, p. 719-724. Discussion by W. F. BERCK 
& J. S. Beaes, v. 72, p. 344-346. 
This paper by the chief engineer of the Marchant Calculating Machine 
Company discusses the engineering questions associated with the mass 
production of calculating machines and similar devices and presents an 














OTHER AIDS TO COMPUTATION 249 


example of the design of such a calculator. The design example includes a 
method for introducing an addend continuously by means of gears into the 
machine. The keyboard selects the proper gear ratio but the actual gear 
shifting is done automatically. There is also a discussion of a carry system 
based on differentials, the use of differentials to position the wheels of a 
register and a continuous method of comparing the remainder in division 
with the divisor in order to control automatic division. F.J.M 


13. B. S. Benson, “Description of an automatic data-analyzing machine,” 

Mech. Engineering, v. 72, 1950, p. 285-288. 

This device has as its input the film record telemetering data from flight 
tests of airplanes and missiles in bar graph form. In each frame the values 
of six out of twenty-eight variables are chosen, measured, and plotted. The 
selection is governed by a counting relay and integrating networks are 
provided for the measurement by photo cells and calibration of the quan- 
tities. The six quantities are distinguished in the graphical recorder by color. 
The time between frames is used as a reference for checking the counter 


relays. F. J. M. 


14. R. L. BisptincHorr, T. H. H. Pian & L. I. Levy, “A mechanical 
analyzer for computing transient stresses in airplane structures,” Jn. 
Appl. Mech., v. 17, 1950, p. 310-314. 

If the point of suspension of a torsion pendulum is rotated in accordance 
with a function ¥(¢), the rotation 6 of the pendulum satisfies the differential 
equation 6 + w?(@ — ¥) = 0 or 6 + w = w(t). Thus a torsion pendulum 
permits an easy realization of this last equation. In the application de- 
scribed here, three such pendulums are set up with the angular output @ 
measured by wire strain gauge elements which when used in bridge circuits 
will yield a.c. voltages proportional to the angle of twist. A linear combi- 
nation of these a.c. voltages may be obtained by using transformers and 
potentiometers. F. J. M. 


15. M. Doporto, ‘Theory and description of a gradient wind computer,” 
Department of Industry and Commerce, Meteorological Service Geo- 
physical Publications, v. 3, No. 1, 1950, 8 p., Stationery Office, Dublin, 
Eire. 

The author devises a graphical computer based on a rather evident 
geometrical construction which solves for the real positive roots of certain 
quadratic equations expressing the meteorological gradient wind velocity. 
Many meteorologists are of the opinion that other approximations (such as 
x = c¢ for large |b|), or nomograms based on a different gradient wind 
equation, give the gradient velocity with sufficient accuracy for most pur- 
poses, considering the approximations inherent in the wind equations used 
and the errors existing in the forecasted values of the parameters involved. 


PAauL STADERMAN 
USN, Special Devices Center 
Port Washington, Long Island, N. Y. 
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16. S. Firer, A. S. MENGEL, G. W. Evans, W. A. Mersman, L. BAvER, 
L. C. MERRILL, W. S. MELaAugN, J. L. Burnsipe, R. N. Brertot, D. L. 
Markusen, C. F. Cook, D. S. TEacuE, R. D. Grtpin, H. H. Goong, 
W. A. WHEATLEY, G. G. DEN BROEDER, R. R. BENNETT, A. S. FULTON, 
H. Exvers, E. VoceL, W. A. McCoo, P. Brock, S. SHERMAN, A. H. 
Miter, H. Metssincer, C. Beck, W. F. RIcHMOND, Symposium on 
REAC techniques, Reeves Instrument Corporation, 215 E. 91 Street, 
New York 28, N. Y., 1951, 148 p. 

Since the war, the use of REAC equipment has become widespread in 
certain fields. A REAC is an electromechanical differential analyzer, manu- 
factured in “packages’’ consisting of seven feedback amplifiers and asso- 
ciated equipment to permit their use as integrators or summers, potentiom- 
eters and four servo units which are used as multipliers, dividers or resolvers. 
The purpose of “Symposium I” was to coordinate the experience of the 
many users of REAC equipment for their mutual benefit. (The symposium 
refers to ‘‘40 installations” many of which have more than one package.) 
The symposium took place on March 15, 16, 1951. 

Session I. The first paper given here was by A. S. Mengel and concerns 
the application of the Rand REAC installation to an unusual type of prob- 
lem in the calculus of variations concerning the motion of a dynamic system 
in which one of the end values is to be maximized. This problem has been 
studied mathematically in previous Rand reports by M. R. HESTENEs who 
showed that its solution depended on the simultaneous solution of the 
equations of motion of the dynamic system and a set of Euler-Lagrange 
equations. The REAC equipment can be effectively utilized both to obtain 
solutions and also to determine the sensitivity of the solutions to various 
parameters which appear. In practical applications, this last is frequently a 
very important question. 

The second paper of this session by G. W. Evans is entitled “Some 
partial differential equations solvable by the REAC.” A single first order 
partial differential equation can be solved by the method of characteristics, 
which is also applicable to certain linear partial differential equations of the 
second order in the hyperbolic and parabolic case. Two examples of these 
last with constant coefficients are given with detailed REAC circuits. 

The third paper by W. A. MERSMAN is concerned with the solution of 
an equation in the form L(x) = M(y) where L and M are linear differential 
operators with constant coefficients and y has a discontinuity at ¢ = 0. 

The fourth paper by L. Bauer & S. FIreEr discusses the realizations of a 
polynomial in a complex variable z by either resolvers which generate sine 
and cosine or by obtaining these functions as the solution of differential 
equations. Either method is applicable to the REAC and such a realization 
of the polynomial can be used to locate the roots. 

The concluding paper of this session is by L. C. MERRILL and is con- 
cerned with the evaluation of multidimensional integrals. These can be 
evaluated as a properly averaged one-dimensional integral on an ergodic 
path through the region of integration. The path is generated by sawtooth 
voltages of incompatible frequencies, and irregular regions of integration in 
the multidimensional space are handled by suitable on-off circuits. The 
range of the integration is extended by a direction reversing method. 
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Session IIT. W. S. MELAHN describes in the first paper of this session 
certain additions to the REAC equipment at Rand including an IBM 
plugboard method of set up, which permits one to retain the set up of an 
individual problem while the machine is used on a different problem, an 
operator’s console, a digital read out and certain stabilized power supplies. 
By additional wiring to the parts of components, the Rand installation has 
attained a much greater flexibility and by additional checking and adjusting 
techniques a greater accuracy than the usual REAC installation. 

The next paper by J. L. BuRNsIDE advocates a number of changes in 
REAC design. One of these is the separation of heat producing non-manipu- 
lated components such as power supplies from the rest, a second is the 
substitution of a switching system for the patch bay set up. He also recom- 
mends the addition of certain equipment like high speed multipliers and 
printing readers for the servo units. 

The third paper by R. N. Bretor & D. L. MARKUSEN is concerned 
with the use of the REAC in the design of the Autopilot for an airplane. 
The design of such a device involves a number of stability considerations 
which are normally considered first in a simplified form by direct calculation 
and then all factors are studied in full generality on the REAC. 

In testing airplane equipment and instruments, it is necessary to con- 
sider their performance under conditions simulating flight. This is usually 
done by mounting these on a rotating mount. C. F. Coox discusses the 
control of one such mount by REAC equipment and finds a region of stable 
control. 

The next paper was by D. S. TEacue & R. D. Gitpim, “Preliminary 
investigation of the suitability of REAC for experimental curve fitting.” 
The objective here is the determination of certain unknown coefficients in a 
differential equation, when a solution is given. 

The concluding paper of this session by H. H. Goopr, W. A. WHEATLEY 
and G. G. DEN BROEDER is concerned with the problem of simulating N 
chance variables with a joint N dimensional distribution with given means 
and given covariance matrix, by analog equipment. This is accomplished 
by means of certain transformations on N variables which are the output 
of N noise generators, whose specifications are given. F.J.M 


Session III. This session consisted of three papers on noise and fre- 
quency analysis and one on the solution of integral equations. 

A. S. FuLton and R. R. BENNETT, “The generation and measurement 
of ultra low frequency random noise.” In this article an electronic noise 
generator suitable for use with an electronic analog computer is described. 
The type of noise produced has a normal amplitude distribution. Its fre- 
quency spectrum is essentially flat at low frequencies with a half power 
point at 25 c.p.s. The mean value is practically zero and is obtained by 
passing the noise through a high pass filter. The frequency spectrum can 
be altered by passing the noise through various types of filters and the 
amplitude distribution can even be changed by suitable combinations of 
linear and non-linear filters. The (residual) mean value of the noise is 
obtained by integrating the voltage for a sufficiently long length of time. 
The mean square value can be determined essentially by half-wave recti- 
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fying the noise and then integrating. The spectrum of a particular noise 
can, of course, be computed from its autocorrelation function. In simulator 
work it is often simpler to pass the noise through various sharp band pass 
filters and measure the transmitted power. 

H. Esters & E. VoGEL, “Simulation of noise in missile homing prob- 
lems.”” (Only an abstract of this paper appears.) It is said to contain a 
description of the noise generator and the shaping of the output to simulate 
the noise as recorded in flight. 

W. A. McCoo, “Frequency analysis by electronic analog methods.” 
In this paper an electronic analog computer technique for evaluating Fourier 
transforms is described. The method is to consider the complex Fourier series 
expansion of period T of an arbitrary function g(é), g(t) = /*. ya(?) 
= >°.. F(wn) exp (jwat). For large values of T it is reasonable to assume 
F(wn) is a good approximation to the Fourier transform of f(t). If, 
in the above equation, F(w,) is replaced by the Fourier transform 
F(w) = f°. f(t) exp (— jwt)dt, there results, on dropping subscripts, 
y(t) = fot f(r) exp (jw(t — r)) dr. From this result, a linear differential equa- 
tion on y(t) (which can be solved by analog techniques) is derived. The 
method also applies to the evaluation of the inverse Fourier transform. 

P. Brock and S. SHERMAN, “Applications of REAC equipment to the solu- 
tion of problems involving integral operators.”’ In this paper practical methods 
(successfully used on REAC equipment) are given for the solution of typical 
problems involving integral operators. The method is by iteration using 
the ‘“‘point-plot method.” For example, the following problem was solved: 
¢(s) = s+ 2 fo! ste(t) dt where go(s) = sand (*) ga(s) = s + 2 fo! sten_s(#) dt. 
¢n(s) is determined from (*) using the previously determined ¢,_; curve as 
a function of t. To solve (*), ga(s) is determined as a point function [¢,(s;) ] 
for a sequence of constant s values [s;]. The REAC solves ¢,(s;, t) = s; 
+ 2 fot sitgn_s(t) dt for each s;. Note that ¢,(s;, 1) = ¢,(s;). This function 
¢n(si) is plotted for each value of s; on an output table. 

Session IV. This session on numerical checks was devoted to a dis- 
cussion by many of the members of the Symposium (under the chairmanship 
of P. Brock) of the general problem of determining the accuracy of REAC 
solutions. No papers were presented. Most of the remarks were of a quali- 
tative nature. 

Session V. This session consisted of five papers devoted to particular 
techniques used on the REAC for solving various special problems. 

A. H. Miter, “A REAC solution.of a linear proton accelerator design 
problem.” The mathematical problem considered in this paper is that of 
solving a number of systems of equations. Each system is composed of one 
linear and one non-linear differential equation. The initial conditions for 
the mth system are determined by the solutions of the (m — 1)st system. 
The problem that arises in solving, these equations on the REAC is that 
seven servos (necessary to generate 1/m, cotangents, tangents, etc.) are 
needed while only four are available. By considering the numerical range 
of the variables involved, certain approximations and other circuit arrange- 
ments can be used. For example, 1/n can be generated by solving the differ- 
ential equation, dy/dn = — y*/100. 

H. MEIssINGER, “Representation of functions of several variables on 
REAC equipment.” One can represent a function of two variables by using 
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a three dimensional cam or a multi-wire pick-off device. Neither of these 
methods is particularly feasible. A method is described where certain func- 
tions can be broken up into sums and products of functions of a single 
variable. For example, it is shown that a certain function 2(x,y) which 
occurred in a recent problem on guided missile trajectory work could -be 
approximated by 2(x, y) = fs(x) — fe(x)fi(y). Approximate analytic expres- 
sions were determined for f,, fe and f; and the result set up on the REAC. 

C. Breck, “A method for solving problems on the REAC by use of 
transfer functions without passive networks.”” Some problems are stated 
directly in terms of the transfer function H(s) of a linear system, that is, 
(*) xo(s) = H(s)x.(s) where s is the Laplace transform variable and H(s) 
is a rational function of s. A method is described whereby (*) can be solved 
without resorting to the original differential equation. The case where the 
system is not originally in the relaxed state (zero initial conditions) is also 
considered. 

W. F. Ricumonp, “REAC techniques.” Various techniques for treating 
continuous linear functions, continuous non-linear functions, discontinuous 
linear functions or linear functions with discontinuous derivatives, discon- 
tinuous non-linear functions or non-linear functions with discontinuous 
derivatives are discussed. Numerous circuit diagrams are included. 

R. R. Bennett, ‘The generation of straight line transfer relationships.” 
The use of relay (sampling, clamping) amplifiers in generating straight line 
segments is described. Such techniques are necessary in the solution of non- 
linear problems on electronic analog computers. 


K. S. MILLER 
New York Univ. 


New York 


17. G. D. McCann & R. H. MACNEAL, “Beam vibration analysis with 
electric-analog computer,”’ Jn. Appl. Mech., v. 17, 1950, p. 13-26. 
Four systems of partial differential equations in two independent vari- 

ables x and ¢ are solved by an analog method in which x is permitted to 

have only a small number of values so that the partial differential equations 
become an ordinary system on the variable ¢, which can be realized by an 
electric analog [MTAC, v. 3, p. 501-513]. The systems treated are linear 
of the second and fourth order on one, two or three unknown functions and 
the coefficients may depend on the variable x. The validity of the results 
also depends on an analogy: results from the analog are compared with 
analytic solutions in certain simple cases. 

F. J. M. 


18. A. D. Moore, “Soap film and sandbed mapper techniques,” Jn. Appl. 
Mech., v. 17, 1950, p. 291-298. 


19. MULTIX Multiplying Sticks. Larchmont, N. Y. $1.00. 

This set consists of ten cardboard strips and two mimeographed pages 
of instruction. The strips are similar to the well known “Napier bones”’ for 
multiplying a number of any size by a one digit number. For each possible 
digit of the multiplicand a strip is given. When these are suitably arranged 
the product can be read off. However in these ‘Multix,” the same strip can 
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be used for either 5 in the multiplicand or 55 or 555 and, consequently, 
only ten strips are needed. The price seems high and many people will find 
multiplication tables more useful, even if less intriguing. 

F. J. M. 


20. D. W. PENDERY, ‘‘604 electronic calculator control panels for multiply, 
divide, and square root on the card-programmed electronic calculator,” 
IBM Tech. Newsletter, No. 1, 1950, p. 4-8. 


W. BELL, ‘The ‘DUZ’ general-purpose control panel for the IBM type 
604, electronic calculator,”” IBM Tech. Newsletter, No. 1, 1950, p. 9-11. 
D. B. MacMrtian & R. H. Stark, “A general-purpose 604 electronic 
calculator control panel for the card-programmed electronic calculator,” 
IBM Tech. Newsletter, No. 1, 1950, p. 12-18. 


B. OLDFIELD, ‘604 electronic calculator diagrams for the calculation of 
sin x, cos x, e*, e~”, sinh x, and cosh x on the card-programmed elec- 
tronic calculator,” IBM Tech. Newsletter, No. 1, 1950, p. 19-25. 

B. G. OLDFIELD, ‘“‘The calculation of variance, covariance, and aver- 
ages,” IBM Tech. Newsletter, No. 2, 1951, p. 2-6. 


B. G. OLDFIELD, ‘‘Numerical integration by Simpson’s rule,” IBM Tech. 
Newsletter, No. 2, 1951, p. 7-12. 

W. D. Bett & Amy McAvory, “Graphical aids for determining shift 
codes,” IBM Tech. Newsletter, No. 2, 195i, p. 13-15. 

D. B. MacMiLian & R. H. Stark, “ ‘Floating decimal’ calculations 
on the IBM card-programmed electronic calculator,” IBM Tech. News- 
letter, No. 2, 1951, p. 16-26. 


The introduction to Technical Newsletter, No. 1, reads as follows: 


“Technical Newsletter Number 1 is the first of a series to be pub- 
lished by the Applied Science Department of the International Business 
Machines Corporation. The purpose of these Newsletters is to allow 
rapid interchange of information between computers concerning com- 
putational methods on IBM punched-card equipment. 

Technical Newsletter Number 1 describes general purpose control 
panel diagrams which have been used successfully. Three of these dia- 
grams are for the IBM Card-Programmed Electronic Calculator. The 
fourth diagram is for the IBM Type 604 Electronic Calculator.” 


In both numbers, the solution to the stated problems is given without 
mathematical discussion in the technical language of the IBM manuals. 
The “shift codes” referred to in the title of the BELL & McAvory paper 
deal with the position of the decimal point. FILM 


21. H. Poritsxy, B. E. Setts & C. E. Danrortn, ‘Graphical, mechanical 
and electrical aids for compressible fluid flow,’ Jn. Appl. Mech., v. 17, 
1950, p. 37-46. 

The flows studied in this paper are two dimensional with no vorticity 
and no viscosity. Consequently, they can be described by two potential 
functions, g and y, by plotting the curves g = C, and wy = C; for equally 
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spaced values of C, and C;. ¢ is the velocity potential and y is the streamline 
function and thus the two sets of curves intersect at right angles. The plot 
divides the given area into curvilinear rectangles, which in the case of 
incompressible flows are equisided. For the compressible case, Bernoulli’s 
equation gives a relation between the sides of these curvilinear rectangles. 
The authors give a historical discussion of graphical methods in this prob- 
lem and describe, with examples, a process of starting with the incom- 
pressible case and arriving at the desired result. For this purpose a number 
of mechanical aids are described. A four legged caliper is described which 
will give the proper size for the other side of the rectangle, when one side 
is given and a cam device for the same purpose is described. To permit a 
ready modification of the plot, the latter is realized by wires and pegs which 
have two holes bored at right angles through their sides. Lord Rayleigh 
has pointed out that the incompressible case can be plotted very effectively 
by using a double mesh of wires. For this purpose each peg has an addi- 
tional pair of perpendicular holes making an angle of 45° with the first pair 
and the double mesh involves both the sides and the diagonals of the smaller 
squares. The authors also discuss the representation of the equations of 
fluid flow by means of a resistance network. Here the resistances are to be 
adjusted until the proper relation between potential and current flow is 
obtained. 
F. J. M. 


22. W. W. Soroka, “Experimental aids in engineering design analysis,” 
Mech. Engineering, v. 71, 1949, p. 831-837. 


This is a survey article with a bibliography, covering certain analogue 
techniques for the solutions of partial and ordinary differential equations 
and algebraic equations. 


23. M. A. Woopsury, ‘“Inverting modified matrices,” Princeton Univ., 
Statistical Research Group, Memorandum Report 42, Princeton, New 
Jersey. 4 hectographed leaves. 

This is a resumé, including the work of SHERMAN and Morrison, of 
formulas for inverting matrices which have been obtained from a matrix 
with a known inverse by bordering or by modifying either an individual 
element or the elements in one row or column or both or in a number of 
rows and columns. ‘‘None of the proofs are difficult but considerable ex- 
ploration was necessary to discover useful forms.” 

F. J. M. 


NOTES 


128. Gauss TO GERLING ON RELAXATION.—In his recent book’ on 
matrices ZURMUHL traces the relaxation method? of solving linear equations 
back to DEDEKIND’s report of Gauss’s lectures.’ It is believed by some 
computers, however, that Gauss’s method is a different one—namely, the 
related method given by SEIDEL.’ In the interest of giving Gauss his proper 
credit as a proponent of relaxation, the following translation of a letter by 
Gauss‘ is offered. Moreover, in this same letter Gauss introduces a useful 
trick,’ now generally forgotten, which seems to improve the convergence 
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of the relaxation process in a wide class of problems. The letter to C. L. 
GERLING dated at Géttingen, 26 Dec. 1823 is as follows.' 


My letter got to the post office too late and was brought back to me. 
I am therefore opening it to add the practical directions on elimination. 
However, this method has many small special advantages which can be 
learned only by use. 

As an example I take your measurements for Orber-Reisig.® 

[The angle measurements communicated by Gerling (taken from a sheet 
found with Gauss’s papers) were as follows, where 1 denotes Berger Warte, 
2 Johannisberg, 3 Taufstein, and 4 Milseburg: 


Repetitions Angle 

13 1.2 = 26°44’ 77423 

28 1.3 = 77 5753.107 

26 1.4 = 136 21 13.481 

50 2.3 = 51 13 46.600 
6 2.4 = 109 37 1.833 
78 3.4 = 58 23 18.161. ] 

I first make 


[the angle toward] 1 = 0, 
whence from 1.3 
3 = 77°57'53"107 


(I prefer this, because 1.3 has more weight than 1.2); then from 








13|1.2/2 = 26°44’ 7423 io eae 
50|2.3|2 = ar? =, AEE 6-088; 
finally from 
26|1.4|4 = 136°21'137481 
6|2.4/4 = 8.529+4 = 136°21'117641. 
78|\3.4|4 = 11.268 
In order to improve the approximation still further, I find, from — 
13}1.2)}1 = — 0°727 
28\1.3;1= 0O 1 = — 07855. 
26|1.4|1 = — 1.840 








Since any common alteration of all directions is allowed when dealing 
only with relative position, I alter ail four by +07855 and set 
1 0° 0’ 07000 + a 
2 26 44 7.551+ 06 
77 57 53.962 + ¢ 
4 = 136 21 12.496 + d. 


In the indirect® process it is very advantageous to ascribe a variation to 
each direction. You can easily convince yourself of this, if you compute the 
same example without this trick,’ so that you lose the great convenience 
of always having as a control the sum of the absolute terms = 0. Now I 
form the four normal equations by the following scheme (in actual applica- 
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tion, and if the terms are more numerous, I separate the positive and 
negative terms) : [The constants are expressed in units of the third decimal 
place. ] 


ab — 1664 ba + 1664 ca + 23940 da — 25610 
ac — 23940 be + 9450 ch — 9450 db + 18672 
ad + 25610 bd — 18672 cd — 29094 de + 29094. 


The normal equations are therefore: 
6+ 67a —13b— 28¢— 26d 
7558 — 13a + 695 — S50c— 6d 
— 14604 — 28a — 50b + 156c — 78d 
+ 22156 — 26a — 6b— 78c + 110d; 
sum = 0. 


coocoo 
+ 


Huu 


In order to eliminate indirectly,* I note that, if 3 of the quantities 
a, b, c, d are set equal to 0, the fourth gets the largest value when d is chosen 
as the fourth.* Naturally every quantity must be determined from its own 
equation, and hence d from the fourth. I therefore set d = — 201 and sub- 
stitute this value. The absolute terms then become: +5232, —6352, +1074, 
+46; the other terms remain the same. 

Now I let 0 take its turn, find b = + 92, substitute, and find the abso- 
lute terms: +4036, —4, —3526, —506. And thus I continue until there 
is nothing more to correct. Of this whole calculation I actually write only 
the following table:* 


d=—201 b=+92 a=—60 c=+12 a=+5 b=-2 a=-1 
+ 6 +5232 +4036 + 16 -320 +415 +41 -—26 


— 7558 -6352 -—- 4 +776 +176 +111 -—27 —-—14 
—14604 +1074 -—3526 -—1846 + 26 -114 -14 +14 
+22156 + 46 — 506 +1054 +118 -— 12 0 +26. 


Insofar as I carry the calculation only to the nearest 2000-th of a second, 
I see that now there is nothing more to correct. I therefore collect 


a=-—60 b=+92 c=+12 d=—201 
+ 5 — 2 
— 1 
—56 +90 +12 —201 
and add the common correction +56, whence 
a= 0 b=+146 c=+68 = —145; 
therefore the values [of the directions ] are 


1 0° 0’ 07000 
2 26 44 7.697 
3 77 57 54.030 
4 136 21 12.351. 


Almost every evening I make a new edition of the tableau, wherever 
there is easy improvement.” Against the monotony of the surveying busi- 
ness, this is always a pleasani entertainment; one can also see immediately 
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whether anything doubtful has crept in, what still remains to be desired, 
etc. I recommend this method to you for imitation. You will hardly ever 
again eliminate directly, at least not when you have more than 2 unknowns. 
The indirect procedure can be done while half asleep, or while thinking 
about other things." 


GEORGE E. ForRSYTHE 
NBSINA 


Univ. of California, Los Angeles 24 


1The preparation of this translation was sponsored (in part) by the Office of Naval 
Research. The author wishes to acknowledge the assistance of Drs. Fritz JoHN and OLGA 
Taussky-Topb. 

2See R. V. SouUTHWELL, Relaxation Methods in Engineering Science. Oxford, 1940, and 
L. Fox, “A short account of relaxation methods,” Quart. Jour. Mech. Appl. Math., v. 1, 
1948, p. 253-280. 

3 Ludwig Seidel, ‘‘Uber ein Verfahren, die Gleichungen, auf welche die Methode der 
kleinsten Quadrate fiihrt, sowie lineare Gleichungen iiberhaupt, durch successive An- 
naherung aufzulésen,” Akad. Wiss., Munich, mat.-nat. Abt. Abhandlungen, v. 11, 1874, 
p. 81-108. Although Seidel’s process is*frequently called the Gauss-Seidel process, I know 
of no place where Gauss mentions it. 

4 The letter is in Gauss’s Werke, v. 9, 1903, p. 278-281. It is probably the letter referred 
to in the footnote, p. 257, of E. T. WHITTAKER and G. Rostnson, The Calculus of Observa- 
ents a London and Glasgow, 1924, as a source of their example of relaxation 

p. 257-8). 

5 Words within brackets are translations of inserts by L. KRUGER, who prepared volume 
9 of Gauss’s Werke. 

6 ‘Indirect elimination’ was Gauss’s term for his iterative process of solving the normal 
equations. It later came to denote any iterative process for solving linear equations. 

7 The symmetric treatment of the unknowns is an essential idea in this and other letters. 
Here Gauss mentions only its advantage as a device which sets up a column-sum check to 
detect errors. In later letters (cf., e.g., Gauss to Gerling, 19 January 1840, Werke, vol. 9, 
p. 250-3) Gauss is convinced that the symmetric treatment of all unknowns yields normal 
equations whose iterative solution converges significantly faster. 

The trick is later described by: Christian Ludwig Gerling (recipient of the letter), Die 
Ausgleichungsrechnung der practischen Geometrie. Hamburg and tha, 1843 (p. 157-8, 
p. 163, p. 386, p. 390); by R. Dedekind, “Gauss in seiner Titeeune iiber die Methode der 
kleinsten Quadrate,”’ Festschrift zur Feter des 150-jihrigen Bestehen der kéniglichen Gesell- 
schaft der Wissenschaften zu Gottingen. Berlin, 1901 (pp. 45-59) and Gesammelte Mathema- 
tische Werke. V. 2, 1931, P- 293-306; and by R. Zurmiihl, Matrizen, Berlin, 1950, p. 280— 
282. (Zurmiihl is wrong, however, in stating that the trick will improve the convergence 
for all badly conditioned systems of equations.) 

For some discussion of when and why the trick may be expected to improve the con- 
vergence of iterative processes for solving linear equations, see George E. Forsythe and 
THEODORE S. Morzkrn, “An extension of Gauss’s transformation for improving the con- 
dition of systems of linear equations,” multilithed typescript at the National Bureau of 
Standards, Los Angeles. 

8 Gauss is here using a method which relaxers recommend: liquidating “that residual 

. which > pope the largest ‘displacement’ ’’ (Fox, op. cit., p. 256). 

® Study of the table shows the algorithm to be precisely the relaxation method described 
by Fox (op. cit., p. 255-6). In comparing Gaiss’s and Fox's presentations, we note that 
Gauss uses the method mentioned in note 8, and that he “liquidates the residuals’’ only 
approximately at each stage, as recommended by Fox (op. cit., p. 257-8) to save unnecessary 
arithmetic. Gauss does not, however, propose ‘‘under-relaxation”’ or ‘“‘over-relaxation,” as 
Fox does. On the other hand, Gauss’s trick mentioned in note 7 is not mentioned by Fox, 
although it is extremely helpful in many common problems. 

10 Gauss wrote, “‘Fast jeden Abend mache ich eine neue Auflage des Tableaus, wo immer 
leicht nachzuhelfen ist.’’ [t is not clear what he meant. 

1! Here one must bear in mind Gauss’s gift for mental arithmetic! 


QUERIES 


49. First USE OF THE TERM HAVERSINE AND THE First TABLE OF 
HAVERSINES.—The editors of the great Oxford New English Dictionary 
(NED) endeavoured to bring together quotations exhibiting the first use 
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in English (so far as known) of each word. In the case of Haversine the 
quotation is from F. G. D. BEprorp (1838-1913), The Sailor’s Pocket Book 

. .» second edition, London, 1875, p. 381. A reference can, however, be 
given to a work published 40 years earlier. In the Dictionary of National 
Biography, v. 10, 1887, J. K. LauGuton’s biography of James INMAN 
(1776-1859), the following sentences occur: “In 1821 appeared his well- 
known book Navigation and Nautical Astronomy for the Use of British 
Seamen, with accompanying tables. In the third edition (1835) he introduced 
a new trigonometrical function, the half-versine, or haversine, the loga- 
rithms of which were added in the tables.’’ Hence here is a new reference, 
to replace that of Bedford, for the editors of NED. The earliest table of 
haversines which I have found is the rather elaborate one in JAMES 
ANDREws’, Astronomical and Nautical Tables, London, 1805 (see MTAC, 
v. 1, p. 422). Can any reader report earlier use of the word haversine, or an 
earlier table of the function? In what Libraries may the 1835 editions of 
Inman’s works (Navigation, and Tables) be found? 

R. C. ARCHIBALD 

Brown University 
Providence, R. I. 


CORRIGENDA 
V. 3, p. 497, 1. —5, for b™ read b™. 
V. 4, p. 124, 1. —2, for 96-97 read 496-497. 
V. 5, p. 25, 854, 1. 7, for 0(.001)25, read 0(.001)16(.01)25. 


’ 
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